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The  objective  of  the  recent  research  effort  was  to  investigate  and  determine  the 
viability  of  utilizing  Locally  Optimal  (L0)  nonlinear  processing  to  mitigate  non- 
Gaussian  interfering  signals  in  a  Direct  Sequence  (DS)  SS  communications  system. 

The  effort  centered  on  the  use  of  memoryless  techniques,  as  well  as  techniques 
employing  memory,  and  performance  comparisons  of  many  receiver  and  nonlinear 
processor  configurations.  The  approach  used  included  the  analysis  and  evaluation 
of  several  implementation  of  the  various  nonlinear  processing  algorithms.  The 
analysis  included  the  study  of  well  known  techniques,  as  well  as  newly  developed 
methods.  Evaluation  was  accomplished  through  the  development  of  software  simulations 
designed  to  test  the  algorithms  in  various  signalling  scenarios.  The  results 
illustrate  the  tradeoffs  of  each  nonlinear  processor  algorithm  for  use  in  a  spread 
spectrum  receiver.  This  knowledge  can  be  used  to  determine  the  most  effective 
processor  for  a  given  interference  scenario.  The  work  presented  in  this  report  is 
directly  in  line  with  the  mission  of  Rome  Laboratory  (RL)  to  provide  secure,  reliable 
communications  to  the  United  States  Air  Force. 
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1.  INTRODUCTION 


Most  modem  military  communication  is  performed  using  Spread  Spectrum  (SS)  systems 
techniques  because  of  their  relative  robustness  in  jamming  environments.  Despite  this  resistance 
to  interfering  signals,  SS  systems  can  undergo  performance  degradation  when  subjected  to  severe 
jamming  [Haze88,  Drap89J.  The  primary  reason  for  this  is  that,  although  SS  techniques  can 
provide  some  immunity  to  jammers,  the  reception  of  the  SS  signal  is  still  accomplished  using 
linear  receivers.  The  linear  receiver  is  only  optimal  in  the  case  of  Gaussian  channel  noise  and 
jamming  signals  are,  by  design,  non-Gaussian. 

It  has  long  been  known  that  the  optimal  and  quasi-optimal  receivers  for  such  scenarios 
are  nonlinear  processors  [Drap85,  Higb88,  Midd66,  Spau85],  These  processors  have,  in  the 
past,  been  difficult  to  implement  in  available  hardware.  Today’s  technology  has  opened  a 
window  of  opportunity  to  solving  the  problem  of  rejecting  interfering  signals  in  such  situations 
because  of  advances  in  high  speed  digital  signal  processors. 

The  objective  of  the  recent  research  effort  was  to  investigate  and  determine  the  viability 
of  utilizing  Locally  Optimal  (LO)  nonlinear  processing  to  mitigate  non-Gaussian  interfering 
signals  in  a  Direct  Sequence  (DS)  SS  communications  system.  The  effort  centered  on  the  use 
of  memoryless  techniques,  as  well  as  techniques  employing  memory,  and  performance 
comparisons  of  many  receiver  and  nonlinear  processor  configurations. 

The  approach  used  included  the  analysis  and  evaluation  of  several  implementations  of  the 
various  nonlinear  processing  algorithms.  The  analysis  included  the  study  of  well  known 
techniques  as  well  as  newly  developed  methods.  Evaluation  was  accomplished  through  the 
development  of  software  simulations  designed  to  test  the  algorithms  in  various  signalling 
scenarios. 
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The  results  illustrate  the  tradeoffs  of  each  nonlinear  processor  algorithm  for  use  in  a 
spread  spectrum  receiver.  This  knowledge  can  be  used  to  determine  the  most  effective  processor 
for  a  given  interference  scenario.  The  work  presented  in  this  report  is  directly  in  line  with  the 
mission  of  Rome  Laboratories  (RL)  to  provide  secure,  reliable  communications  to  the  United 
States  Air  Force. 


1.1  Problem  Statement  and  Formulation 

The  design  and  development  of  a  digital  communications  system  which  exhibits  a  low 
probability  of  error  (PJ  in  a  variety  of  jamming  scenarios  is  of  paramount  importance  in  military 
applications.  To  attain  this  goal,  it  is  necessary  to  design  a  receiver  which  will  mitigate  the 
effects  of  the  interference  and  thereby  facilitate  the  correct  demodulation  of  the  transmitted 
signal.  For  many  practical  applications  in  which  the  interference  is  thermal  (Gaussian1)  noise, 
the  optimum  receiver  is  a  linear  processor  which  incorporates  the  use  of  a  matched  filter  for 
signal  recovery.  The  matched  filter  impulse  response  is: 

h{t)=ks(T-t)  (1_1) 

where  h(t)  is  the  matched  filter  impulse  response,  s(t)  is  the  real  transmitted  signal  waveform, 
T  is  the  sampling  time  of  the  matched  filter  (equal  to  the  period  of  the  transmitted  waveform), 
and  k  is  an  arbitrary  constant. 

Alternately,  a  time  correlator  realization  of  the  matched  filter  can  be  used  for  signal 
recovery.  The  time  correlator  produces  an  output  exactly  equal  to  that  of  the  matched  filter  at 


1  Thermal  noise  can  be  modelled  using  the  Gaussian  Probability  Density  Function  (PDF)  given  by: 
j  (n-4)1 

f^n)  =  —  e  ,  where  N  is  the  Gaussian  random  variable,  p  is  the  mean,  and  a2  is  the  variance. 
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the  sample  time  T,  and  can  be  written  as: 


T 

z=fr{t)s(t)dt 


0 


(1-2) 


where  z  is  the  output  of  the  correlator  sampled  at  time  T,  r(t)  is  the  received  signal,  and  s(t) 
is  the  transmitted  signal  waveform.  In  most  digital  communications  systems,  the  received  signal 
is  sampled  and  discrete-time  processing  is  used.  Therefore,  Eq.  (1-2)  becomes 


l 


.-E 


m*l,...,Af 


(1-3) 


where  4,  is  the  correlator  output  corresponding  to  the  m*  possible  transmitted  signal,  rk  is  the 
Id1'  sample  of  the  received  signal,  is  the  k*  sample  of  the  rrf1  possible  transmitted  signal,  N 
is  the  number  of  samples  (elements)  per  signal  vector,  and  M  is  the  total  number  of  possible 
signals.  As  may  be  seen,  Eq.  (1-3)  is  a  linear  function  of  the  received  samples,  and  hence  is 
termed  a  linear  receiver.  Note  that  when  the  communications  system  employs  two-dimensional 
signalling,  the  operations  of  Eqs.  (1-1)  to  (1-3)  must  be  performed  on  both  the  In-phase  (I)  and 
Quadrature  (Q)  channels. 


In  a  jamming  environment,  however,  a  digital  communications  system  that  employs  a 
linear  receiver  experiences  a  performance  degradation  due  to  the  highly  non-Gaussian  nature  of 
the  jammer.  Therefore,  communications  systems  subjected  to  large  jamming  signals  require 
other  (nonlinear)  processing  methods  for  signal  recovery.  Furthermore,  the  receiver  must  be 
adaptive  since  in  hostile  situations  the  jamming  statistics  will  often  be  nonstationary. 

The  nonlinearity  needed  can  be  derived  from  the  likelihood  function  of  statistical  decision 
theory  [Midd66,Spau85].  In  the  classical  approach  two  methods  are  discussed:  they  are  the 
Globally  Optimal  (GO)  and  the  Locally  Optimal  (LO)  demodulators.  The  GO  demodulator 
provides  more  robust  performance  whereas  the  LO  demodulator  is  much  simpler  to  implement 
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[Haze88].  The  former  is  obtained  directly  from  the  likelihood  function  whereas  the  latter  is  an 
alternate  form  obtained  from  a  Taylor  series  approximation  and  yields  near-optimal  performance 
for  large  Jammer-to-Signal  ( J/S)  ratios.  The  likelihood  function  is  generated  from  the  PDF  of 
the  incoming  signal. 


The  memoryless  two-dimensional  GO  demodulator  decision  rule  is  formulated  in  the 


following  manner. 

For  a'received  signal  vector  pair  ( rn  rQ) : 


Choose  (s,  ,sQ  )  which  maximizes 


,v 

*•1  L 


(1-4) 


where  ,  sQ )  is  the  m *  of  M  possible  transmitted  signal  vector  pairs  with  components  (s^Sg)  , 
(r,,rn)  is  the  Jfc*  sample  pair  of  the  received  signal  vector  pair,  N,  and  NQ  are  the  I  and  Q 
random  noise  vectors,  and  fNN(  • )  is  the  noise  PDF. 


The  memory  less  two-dimensional  LO  demodulator  is  obtained  from  the  GO  demodulator 
through  a  Taylor  Series  expansion  around  the  signalling  point.  The  resulting  decision  rule  is: 
Choose  (Ij ,  sQ )  which  maximizes: 


.V 

/*,=£  fe-  ^(r,)cos0,+5o^(rt)sind 

t»i 


(1-5) 


where  is  the  envelope  of  the  received  signal  sample  pair,  0t=tan 


-i 


c. 


is  the 


phase,  fR(r A)  is  the  received  envelope  PDF,  and 
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£('■*)*- 


~fR(rk) 

<*r*  \  1 

/«('*)  *  rk 


is  the  LO  Memory  less  Nonlinear  Transform  (MNT). 


(i-tf) 


Both  the  GO  and  LO  demodulators  require  knowledge  of  the  received  signal  PDF.  In 
addition,  the  LO  demodulator  uses  the  derivative  of  the  received  signal  PDF.  Thus,  both 
receivers  need  a  mechanism  to  derive  (an  estimate  of)  the  PDF  from  the  input.  Historically, 
experimental  PDFs  have  been  estimated  using  histograms  with  bins  of  equal  width.  One 
modification,  called  the  Equi-probable  Bin  Histogram  (EBH)  method  [Haze88],  uses  bins 
containing  equal  numbers  of  samples. 

One  major  drawback  inherent  to  the  histogram  methods  of  estimating  the  received  signal 
PDF  is  that  they  produce  discrete  approximations  of  the  actual  continuous  PDF.  The  error 
associated  with  these  approximations  is  further  amplified  by  the  differentiation  process  required 
to  produce  the  LO  MNT.  Thus,  a  means  of  constructing  a  continuous  estimate  of  the  received 
signal  PDF  would  be  beneficial.  One  such  method  which  provides  a  piecewise  continuous 
approximation  is  the  M-Interval  Polynomial  Approximation  (MIPA)  [Haze88].  However,  this 
approximation  still  has  discontinuities  at  the  breakpoints  where  adjacent  approximating 
polynomials  meet.  The  Fourier  Series  Approximation  (FSA)  approach,  on  the  other  hand, 
provides  an  estimate  of  the  PDF,  using  the  Fourier  series  expansion,  which  is  continuous 
throughout  its  entire  domain. 

The  derivation  and  performance  evaluation  of  the  EBH  and  the  FSA  memoryless  LO 
processor  implementations  have  comprised  a  major  portion  of  the  recent  research  effort.  The 
simulation  of  both  memoryless  methods  in  a  standard  communications  system  and  also  a  spread 
spectrum  system  provided  valuable  insight  into  the  LO  processor’s  performance  characteristics. 
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Also,  the  performance  of  FSA  and  EBH  LO  processing  in  a  spread  spectrum  system  was 
compared  to  the  performance  associated  with  standard  histogram  and  MIPA  processing.  Finally, 
LO  processing  with  memory  was  examined,  with  emphasis  on  the  methodology  incorporated  in 
the  Robust  Digital  Adaptive  Transceiver  [Drap89]. 


1.2  Development  of  the  Memoryless  Locally  Optimal  (LO)  Maximum  Likelihood  (ML) 
Function  for  Two-Dimensional  Signalling 

A  brief  derivation  of  the  memoryless  GO  and  LO  detectors  is  presented  in  this 
subsection. 

Let  the  transmitter  be  a  two-dimensional  M-ary  signal  source  represented  by  a  random 
vector  pair  with  an  I  component  S,  taking  on  a  value  s,  ,  and  a  Q  component  50  taking  on  a 

value  sQ  .  Here  (s,  ,sQ)  represents  the  m*  of  M  possible  signal  pairs.  For  the  received  signal, 

the  I  component  is  R,  with  value  r,  and  the  Q  component  is  ~RQ  with  value  ~rQ.  Finally  the 

interference  (thermal  noise  and  jammer)  is  characterized  by  an  I  component  N,  having  a  valuer^ 

and  a  Q  component  NQ  having  a  value  nQ.  The  relationship  between  the  transmitted  and 
received  signals  is  given  by: 

*,-P,  .  (1'7) 

The  Maximum  Likelihood  (ML)  detection  algorithm  at  the  receiver  is  implemented  as: 
Choose  the  possible  signal  pair,  (s,  ,sg),  which  maximizes: 
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i 


««*/**(  (1*8) 

which  is  the  PDF  of  the  received  signal  conditioned  on  the  m*  transmitted  signal  pair.  The 
variable  zm  is  called  the  ML  decision  statistic.  Using  the  relationship  defined  in  Eq.  (1-7)  and 
assuming  statistical  independence  of  the  signal  and  interference,  Eq.  (1-8)  becomes: 

ri~  Si.,ra~  sqJ  (1-9) 


where  is  the  joint  interference  PDF.  Next,  if  the  elements  of  both  the  I  and  the  Q 

received  vectors  are  assumed  to  be  Independent  and  Identically  Distributed  (IID)2,  Eq.  (1-9) 
becomes 


jfc*l 


(1-10) 


where  (r,,rQ)  is  the  k*  received  signal  sample  pair,  and  ( s ,  ,sQ  )  is  the  k*  sample  pair  of  the 

m*  possible  source  signal.  Since  the  logarithm  is  a  monotonically  increasing  function,  the 
natural  logarithm  of  Eq.  (1-10)  is  used  to  simplify  calculations  [Mels78].  This  results  in  the  GO 
ML  decision  statistic  of  Eq.(l-4),  repeated  here  for  convenience: 

Choose  (^,s0)  which  maximizes: 


ln[^Wr/,  ~SL’ru. 


(Ml) 


:  The  resulting  GO  and  LO  processors  are  memoryless  due  to  the  IID  assumption. 
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For  a  large  J/S  ratio3,  a  first-order  two-dimensional  Taylor  Series  expansion  around  the 
received  signal  point  can  be  used  to  approximate  the  interference  PDF.  The  expansion  is  valid 
under  the  assumption  that  the  deviation  of  the  interference  from  the  received  signal  point  is 
minimal  for  a  large  J/S.  Also,  in  this  case  the  interference  PDF  may  be  approximated  by  the 
received  signal  PDF,  fRR( • ) .  These  simplifications  result  in  the  LO  ML  decision  statistic 

[Illi91]: 


Choose  ( Sj ,  sQ )  which  maximizes: 


N  < 


'„-E 


*•1 


37 -A*Sri>’rQ)  Jpf*fiSri,'ra) 

_Ji _ g| _ 

“  fn,R Sri,’rQ) 
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However,  determination  of  the  received  signal  joint  PDF  is  a  difficult  process,  and  determination 
of  the  partial  derivatives  even  more  so.  To  reduce  the  complexity  of  the  likelihood  function, 
bivariate  radial  symmetry  of  the  received  signal  PDF  is  assumed.  Under  this  assumption,  the 
received  signal  PDF  is  given  by: 


Ur) 

2irr ' 

0, 


O<0<2t 

elsewhere 


(M3) 


where  fR(r)  is  the  received  envelope  PDF,  r=\Jrf+r2Q  is  the  received  magnitude,  and 
0=tan"‘(rQ/r,)  is  the  received  phase  angle.  The  radial  symmetry  assumption  is  valid  for 


3  A  typical  jammer  is  usually  at  least  20  dB  greater  than  the  transmitted  signal.  Therefore,  for  a  signal 
environment  where  jammers  are  present,  a  large  J/S  assumption  is  valid. 
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interference  sources  of  interest  as  they  will  have  random  phase  angles.  Even  a  constant- 
frequency  waveform  will  have  a  vector  that  rotates  at  a  uniform  rate  and  is  therefore  equally 
likely  at  any  angle  [Higb88].  With  this  assumption  the  LO  ML  decision  statistic  of  Eq.  (1-12) 
reduces  to  Eqs.  (1-5)  and  (1-6),  repeated  here  for  convenience: 

Choose  (s,  ,sQ)  which  maximizes: 

{sLg(rk)cosdk+sQj(rk)sinOk}  (1'14) 

where 


is  the  LO  MNT. 


(1-15) 
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2. 


MEMORYLESS  LOCALLY  OPTIMAL  PROCESSOR  ALGORITHMS 


The  memory  less  LO  receiver  derived  in  Section  1.2  requires  the  PDF  of  the  received 
signal  to  generate  the  LO  MNT  used  in  detection  of  the  transmitted  signal.  However,  in  practice 
a  priori  knowledge  concerning  the  received  signal  PDF  is  not  available.  Instead,  an 
approximation  must  be  made  based  on  the  information,  usually  the  received  signal  samples, 
available  at  the  receiver.  This  section  presents  two  methods  of  estimating  the  required  PDF. 

The  Equiprobable  Bin  Histogram  (EBH)  method  is  an  extension  of  the  traditional 
histogram  method  [Krey88]  of  estimating  PDFs.  In  this  method,  the  histogram  bins  are  chosen 
to  have  equal  probability  rather  than  equal  width  as  is  the  usual  case.  Thus,  the  width  of  the 
bins  is  allowed  to  vary,  while  the  height  of  each  bin  is  constrained  so  that  the  area  under  the 
PDF  approximation  is  equal  to  one.  This  provides  an  improved  resolution  over  the  traditional 
histogram  method  in  regions  where  there  are  a  large  number  of  samples. 

The  next  method  discussed  is  the  Fourier  Series  Approximation  (FSA)  which  utilizes  a 
Fourier  series  to  estimate  the  received  signal  PDF.  This  algorithm,  in  contrast  to  the  EBH 
method,  produces  an  approximation  of  the  PDF  which  is  continuous  and  has  continuous 
derivatives  throughout  its  entire  domain.  This  is  extremely  useful  since  the  LO  MNT  requires 
the  derivative  of  the  received  signal  PDF,  and  the  differentiation  process  tends  to  increase  the 
error  associated  with  the  discrete  approximations  generated  by  the  histogram  methods. 
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2.1  Equiprobable  Bin  Histogram  (EBH) 

A  standard  histogram  has  bins  of  equal  width,  and  provides  equal  resolution  over  the 
entire  range  of  interest.  Intuitively  it  seems  better  if  more  resolution  could  be  given  to  the 
intervals  of  data  corresponding  to  higher  probabilities.  One  method  of  achieving  this  is  to 
constrain  the  bins  to  have  equal  probability  rather  than  equal  width.  The  probability  of  a 
random  sample  falling  in  a  given  bin  is  the  area  of  the  bin,  so  the  bins  must  have  equal  area  to 
be  equiprobable.  For  a  given  data  sequence  this  may  be  accomplished  by  assigning  an  equal 
number  of  samples  to  each  bin.  A  histogram  with  these  constraints  is  known  as  an  Equiprobable 
Bin  Histogram  (EBH),  and  unlike  a  standard  histogram,  both  the  height  and  the  width  of  the 
EBH  bins  vary. 

An  EBH  is  constructed  by  sorting  the  data  in  ascending  order  and  assigning  an  equal 
number  of  data  points  to  each  bin.  Figure  (2-1)  is  a  five  bin  (B~ 5)  EBH  for  an  arbitrary  data 
sequence.  The  data  point  with  the  largest  value  in  the  i*  bin  is  called  the  breakpoint  and  is 
denoted  X(.  The  zeroth  breakpoint,  X0,  is  the  minimum  value  of  all  data  points.  The  height  of 
each  bin,  Y„  is  determined  from  the  constraints  described  above:  the  width  times  the  height  of 
each  bin  must  be  equal  and  the  total  area  of  all  bins  must  equal  one.  Thus,  Yt»(Xi-Xt.l)/B. 

Five  Bin  EBH 


Amo i . tuae 

-20  CO  20  -JO  60  80 


am,  •  ;  a ;  »  -  :  * 


Figure  (2-1) 
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Several  examples  demonstrate  the  differences  between  a  standard  histogram,  an  EBH, 
and  the  ideal  PDF.  A  6  dB  Binary  Phase  Shift  Keying  (BPSK)  signal  added  to  0  dB  Gaussian 
noise  is  shown  in  Fig.  (2-2).  Figure  (2-3)  shows  the  corresponding  10  bin  histogram  and  10  bin 
EBH.  As  may  be  seen  in  the  figure,  both  histogram  methods  provide  an  approximation  that  is 
reasonably  close  to  the  ideal  PDF,  represented  by  the  smooth  non-boldface  line.  The  same 
analysis  for  a  26  dB  BPSK  signal  added  to  0  dB  Gaussian  noise  is  shown  in  Fig.  (2-4).  In  this 
case  there  are  two  regions  of  high  probability,  namely  at  amplitudes  of  -20  and  20.  It  may  be 
seen  that  the  standard  histogram  gives  equal  resolution  to  the  entire  range  of  the  PDF,  while  the 
EBH  focuses  on  the  areas  of  high  probability  and  provides  a  representation  that  is  closer  to  the 
ideal  PDF. 


First  lOOO  Samples  of  BPSK  Data 


O  200  400  600  600  lOOO 

Sample 

Figure  (2-2) 


It  should  be  noted  that  the  magnitude  PDFs  of  the  jamming  scenarios  in  this  report  have 
a  single  region  of  high  probability,  and  the  performance  of  the  EBH  offers  little  or  no 
improvement  over  the  standard  histogram.  However,  for  future  jamming  scenarios  the  EBH 
may  be  of  greater  value. 
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Probtb 


2.2  The  Fourier  Series  Approximation  (FSA)  Algorithm 

A  novel  method  of  generating  a  continuous  LO  MNT  using  the  Fourier  series  has  been 
derived  and  simulated.  The  Fourier  Series  Approximation  (FSA)  is  a  method  of  generating  a 
continuous  estimation  of  the  PDF  of  a  random  variable.  This  method  is  extremely  useful  for 
LO  processing  since  the  LO  MNT  requires  the  derivative  of  the  received  signal  PDF,  and  the 
differentiation  process  amplifies  the  error  associated  with  discrete  histogram  approximation 
methods.  This  section  presents  the  derivation  of  the  FSA  algorithm,  a  method  of  implementation 
based  on  a  histogram  approach,  and  a  direct  parametric  implementation  of  the  FSA  algorithm. 


2.2.1  Derivation  of  the  FSA  Algorithm 

Given  the  received  signal  magnitude,  r,  the  envelope  PDF,  fR(r),  is  approximated  by  the 
following  expression: 

f*(r) = ^  ancos(n(j0r)  +£  bnsin(no>0r)  (2-1) 

n»  I  n»l 


2  -jj~ 

where  o>0  =  - - —  =  — ,  is  the  maximum  value  of  the  received  signal  set  being  processed,  rm 


max  mm 


is  the  minimum  value  of  the  received  signal  set,  and  p  is  the  order  of  the  FSA.  The  coefficients 
{a„}  and  {bn}  are  determined  by  minimizing  the  following  error  criterion: 
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(2-2) 


r 


C4(r)-^-£  a„cos(mv)-£  Mn(/I<V)]2<*’ 

•  *«i  *»i 


Differentiating  Eq.  (2-2)  with  respect  to  a#  and  equating  to  zero  yields  the  following  expression: 


f  4  +  E  a.cosCncj/)  +  ^sin(/iWor)  dr  »  f  /*(r) dr 

)  I  ^  *«i  «-t  / 


By  noting  that  the  integration  of  the  sine  and  cosine  functions  over  an  integer  multiple  of  their 

r _ 

period  is  equal  to  zero  and  that  [  fR{r)dr-  l ,  Eq.  (2-3)  becomes 


(tux  tnm 


(2-4) 


To  find  the  coefficients  {a,},  differentiation  is  performed  with  respect  to  ak  and  the  result  is 
equated  to  zero.  This  gives  the  following  equation: 


f  cos(kw0r)  fR{r)  -  J!  -  53  a„cos(nw0r)  -  53  &„sin(itoy)  Jr = 0 

/  *  n*l  n»l 


After  some  manipulation,  the  expression  for  ak  is  found  to  be 

rjt 

ak= - 1 — f  cos(ku0ryR(r)dr,  k=\,...,p  (2’& 

^"max  ^"min  / 

Similarly,  the  coefficients  {&„}  are  found  by  differentiating  Eq.  (2-2)  with  respect  to  bk  and 
equating  to  zero,  resulting  in  the  expression: 
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j  sin(*<j0r) 


fR(r)  - ^  ’E  fl«cos (no/)  -£  d.sin^u/) 

*  *•! 


/i«l 


dr=Q 


Again,  after  some  manipulation,  the  expression  for  bk  is  found  to  be 


'« 


bk= - f  sin (kw0r)fR(r)dr,  k=l,...,p 

^aux  ^aim 


(2-7) 


(2-8) 


Finally,  the  LO  MNT  function  is  given  by: 


g(r)=- 


>>  1 
m  r 


(2-9) 


and  its  FSA  derived  estimate,  g(r ) ,  is 


y  y 

£  ^nwoSinf/iu^)-^  bnu0cos(nu</) 

n«l  n»l  .  1 

_ . _  +  _ 

a  p  p  r 

4+Y,  a„cos(/tay)+E  6„sin(nco0r) 

*•  n»l  n»l 


(2-10) 


Equations  (2-1),  (2-4),  (2-6),  (2-8),  and  (2-10)  specify  the  FSA  algorithm  completely. 


2.2.2  The  Histogram-based  Implementation  of  the  FSA  Algorithm 

One  method  of  implementing  the  FSA  algorithm  involves  using  the  FSA  to  approximate 
a  histogram  estimate  of  the  PDF.  The  coefficients  {a„}  and  {bn}  are  approximated  by  evaluating 
Eqs.  (2-6)  and  (2-8)  using  rectangular  numerical  integration  [Krey88],  with  fR(r)  replaced  by  a 
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histogram  approximation  of  the  received  signal  envelope  PDF.  The  resulting  approximations 
for  the  FSA  coefficients  are: 
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£  y,cos^n 

i*0 


(2-11) 


E  »-12> 

max  'am  ‘*0  'max  'mm 

where  yi  is  the  value  of  the  t*  histogram  bin  of  the  PDF,  x,  is  the  received  signal  value 

corresponding  to  the  i*  histogram  bin,  B  is  the  number  of  histogram  bins,  and  A  is  the  bin 
width. 


The  LO  MNT  is  then  estimated  using  the  following  algorithm: 

1.  Let  N,  =  N2  =  D2  =  0  and  let  D,  =  -y  :  (INITIALIZATION) 

2.  For  n  =  1  to  p  :  (ITERATION) 

{ 

Nx  =/V,  *o*w0sin(/7oy) 

N2=N2~  bnnu0cos(nu0r) 

D,  =D,  +  dncos(/iw0r) 

D2=D2+bnsin(n<o0r) 

} 

N.+N,  1 

3.  g(r)  =  -J—J.+- i  :  (COMPUTATION) 
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2.2.3  Parametric  Implementation  of  the  FSA  Algorithm 

The  expectation  of  a  function  h (X)  of  a  random  variable  X  having  a  PDF  fx(x )  is  defined 
as 

OB 

E{h(X)}=ih(x)f^x)dx  (2-13) 

Using  this  definition,  Eqs.  (2-6)  and  (2-8)  may  be  rewritten  as 

a  * — ? — EJcos^Wq/?)},  (2-14) 

r-r 

max  mm 

and 

bn  = - - - Efsin^Oo/?)},  n=l,...,p  (2-15) 

^mu 


Equations  (2-14)  and  (2-15)  can  be  approximated  using  a  sample  average  of  the  form 


a.-— 1—  4Ec°s(W.) 

rwx  rmm  ™  *•' 


(2-16) 


1 

»"(/!«/*) 

rv  *-i 


(2-17) 


where  rk  is  the  k*  sample  of  the  random  variable  R,  and  N  is  the  total  number  of  available 
samples. 


The  parametric  implementation  of  the  FSA  has  an  advantage  over  the  histogram-based 
method  in  that  the  parametric  method  does  not  require  the  intermediate  stage  of  approximating 


the  PDF  with  a  histogram.  However,  the  parametric  form  has  the  disadvantage  of  added 
complexity  since  the  sample  averages  of  Eqs.  (2*16)  and  (2-17)  require  calculation  of  2Np 
trigonometric  functions,  where  p  and  N  can  both  be  large. 
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3. 


SIMULATION  OF  THE  MEMORYLESS  LOCALLY  OPTIMAL  PROCESSOR 
ALGORITHMS 


Software  simulation  was  the  method  chosen  to  characterize  the  performance  of  each 
memoryless  LO  receiver  implementation.  For  most  types  of  jammers,  determining  the 
theoretical  performance  of  the  LO  receiver  is  a  mathematically  intractable  problem.  Simulation 
provides  a  means  of  determining  the  types  of  performance  gains  incurred  in  a  variety  of  jamming 
scenarios. 

Two  types  of  communications  systems  were  simulated  using  the  Signal  Processing 
Worksystem™  (SPW)  by  Comdisco  Systems,  Inc.  The  First  was  a  conventional  Quadrature 
Phase  Shift  Keying  (QPSK)  communications  system.  This  simulation  was  used  to  isolate  the 
performance  gain  associated  with  the  LO  processor.  The  second  simulation  configuration  was 
a  QPSK  DSSS  system  which  was  used  to  determine  the  performance  characteristics  of  the  LO 
processor  in  a  spreading  environment.  By  using  both  simulations,  it  is  possible  to  determine 
scenarios  in  which  the  LO  processor  improves  overall  system  performance. 

This  section  provides  an  introduction  and  overview  of  the  SPW  simulation  package.  The 
conventional  QPSK  and  QPSK  DSSS  simulations  are  also  presented. 
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3.1  Introduction  to  the  SPW  Environment 


The  Signal  Processing  Worksystem™  (SPW)  from  Comdisco  Systems,  Inc.  is  the 
simulation  platform  of  choice.  It  consists  of  several  modules  which  aid  in  the  different  design 
phases  of  communication  systems.  These  include:  the  Block  Diagram  Editor  (BDE),  the  Signal 
Display  Editor  (SDE),  the  Filter  Design  System  (FDS),  and  the  Code  Generation  System  (CGS) 
(Comd91). 

The  BDE  is  a  graphical,  block  oriented  environment  in  which  communication  systems 
can  be  readily  implemented  and  simulated.  Each  system  component  is  represented  as  a  graphical 
block.  The  blocks  provided  in  SPW  are  very  general,  allowing  use  in  a  wide  variety  of 
applications.  Over  150  library  blocks  are  provided  in  the  BDE,  including: 

•  modulators  /  demodulators  such  as  Binary  Phase  Shift  Keying  (BPSK),  Quadrature  Phase 
Shift  Keying  (QPSK),  and  Minimum  Shift  Keying  (MSK) 

•  encoders  /  decoders  such  as  Hamming,  Reed-Solomon,  convolutional 

•  mathematical  functions,  including  simple  arithmetic,  trigonometric,  logarithmic, 
integration  and  differentiation 

•  filters,  including  Butterworth,  Chebyshev,  Elliptic,  Bessel,  and  user  defined  filters. 

•  vector  functions,  such  as  inner  products,  serial-to-parallel  and  parallel-to-serial 
conversions 

•  signal  sources  /  sinks,  including  binary  sources,  white  noise  generators,  pseudo-random 
number  generators,  and  user  defined  sources. 

In  addition  to  these  blocks,  the  user  may  define  new  blocks,  either  through  a  hierarchical 
combination  of  existing  blocks  or  by  designing  custom  coded  blocks  in  Fortran  or  the  C 
Programming  Language  (C). 

The  SDE  can  be  used  to  generate  signals  for  use  in  the  BDE,  or  to  examine  the 
simulation  results  from  a  system  that  has  been  successfully  modeled  and  simulated  in  the  BDE. 
The  major  strength  of  the  SDE  lies  in  its  signal  analysis  capabilities.  In  particular,  the  following 


3-2 


functions  are  available:  Fast  Fourier  Transforms  (FFTs),  histograms,  autocorrelation  and 
cross-correlation  analysis,  and  spectral  density  evaluation.  Also,  signals  may  be  directly 
modified  by  a  host  of  operations.  For  example,  addition,  scaling,  inversion,  differentiation, 
integration,  and  trigonometric  functions  may  be  performed  on  a  signal  or  group  of  signals.  In 
addition,  several  nonlinear  filtering  operations  are  provided,  including  quantizing,  clipping,  and 
limiting. 

The  FDS  allows  the  designer  to  model  and  study  many  different  types  of  digital  filters 
for  use  in  communication  systems.  Any  Finite  Impulse  Response  (FIR)  or  Infinite  Impulse 
Response  (HR)  filter  can  be  designed.  Butterworth,  Chebyshev,  Elliptic,  Equiripple,  and  Bessel 
filters  may  be  easily  implemented  with  the  aid  of  pre-defined  "methods".  In  addition,  the  user 
can  define  methods  for  any  other  type  of  filter.  Once  a  filter  is  modeled  with  the  FDS  it  can 
be  used  in  the  BDE  or  SDE. 

The  CGS  is  able  to  convert  the  BDE  simulations  to  C  code.  This  code  can  then  be 
downloaded  to  any  of  the  digital  signal  processing  boards  or  computers  supported  by  SPW. 


3.1.1  Classification  of  Blocks 

The  blocks  in  the  BDE  may  be  classified  into  three  categories:  SPW  primitive  and 
hierarchical  library  blocks,  user  custom  coded  primitive  blocks,  and  user  hierarchical  blocks. 
Library  blocks  are  the  blocks  supplied  as  a  part  of  the  SPW  environment.  Some  library  blocks 
are  primitive,  in  that  their  specification  is  directly  linked  to  the  underlying  software  code  (e.g 
C).  Others  are  hierarchical,  meaning  that  they  are  composed  of  other  SPW  library  blocks  which 
may  be  either  primitive  or  hierarchical.  A  hierarchical  block  is  specified  by  means  of  a 
"detail"1  which  displays  the  interconnection  of  its  constituent  blocks.  User  custom  coded 


1  Words  enclosed  in  quotes  ("")  refer  to  SPW  terminology. 
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primitive  blocks  are  written  by  the  user  in  C  or  Fortran  and  linked  to  a  graphical  symbol.  User 
hierarchical  blocks  are  made  up  of  library,  custom  coded,  or  other  user  hierarchical  blocks.  As 
may  be  seen,  the  capacity  to  create  multiple  levels  of  hierarchy  provides  an  extremely  powerful 
and  flexible  simulation  environment. 

Figure  (3- la)  shows  the  symbol  for  a  user  custom  coded  primitive  bubble  sort  block,  and 
Fig.  (3- lb)  shows  its  detail.  The  detail  of  this  block  is  a  list  of  parameters  that  are  passed  to 
the  C  or  Fortran  code  when  the  system  is  simulated.  Once  the  block  is  complete,  the  SPW  user 
need  not  be  concerned  with  the  actual  C  or  Fortran  code. 

Figure  (3-2a)  shows  the  symbol  for  a  user  hierarchical  MNT  block,  the  detail  of  which 
is  displayed  in  Fig.  (3-2b).  The  inputs  and  outputs  of  the  symbol  in  Fig.  (3-2a)  correspond  to 
"ports"  in  the  detail  in  Fig.  (3-2b).  The  ports  must  have  the  same  names  as  the  inputs  and 
outputs,  but  need  not  be  in  the  same  position.  Note  that  in  addition  to  the  parameter  list,  other 
blocks  are  included  in  the  detail.  This  is  the  characteristic  difference  between  hierarchical  and 
primitive  blocks. 

The  concept  of  hierarchy  is  very  useful  in  system  design,  and  is  an  integral  part  of  the 
BDE.  On  the  top  level  of  the  hierarchy  is  the  system  to  be  simulated.  The  system  consists  of 
several  blocks.  These  blocks  may  contain  sub-systems  (hierarchical  blocks),  which  are  a  lower 
level  of  hierarchy.  At  the  lowest  level  of  the  hierarchy  are  blocks  that  point  directly  to  C  or 
Fortran  code,  i.e.  primitive  blocks.  One  useful  aspect  of  hierarchy  is  that  a  system  can  be  easily 
transformed  into  a  hierarchical  block.  This  block  may  then  be  used  as  a  component  of  a  larger 
system,  and  the  process  continued  until  the  overall  system  is  completed. 

Experience  has  shown  that  judicious  use  of  all  three  types  of  blocks  provides  the  most 
efficient  simulations,  both  in  terms  of  system  complexity  and  simulation  execution  time.  By 
taking  advantage  of  the  hierarchical  nature  of  SPW,  systems  can  be  modularized,  greatly 
simplifying  their  analysis. 
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(a) 


Figure  (3-1) 


Poly  UNT  Block 


3.1.2  Exported  Parameters 


A  very  useful  feature  of  the  BDE  is  the  exportation  of  parameters.  This  enables  the 
designer  to  set  the  value  of  a  parameter  equal  to  a  function  of  other  parameters.  Often,  several 
different  blocks  have  a  common  parameter,  and  this  allows  all  common  parameters  to  be 
exported  to  a  single  parameter  on  the  top  level  of  hierarchy.  Rather  than  "pushing"  into  each 
block  and  changing  each  parameter,  the  designer  need  only  change  the  parameter  on  the  top 
level  and  all  of  the  exported  parameters  will  assume  that  value. 

3.1.3  On-Line  Documentation 

On-line  help  may  be  obtained  in  the  BDE  by  selecting  HELP  on  the  EZ-Menu  [Comd91] 
and  then  pointing  to  the  desired  object.  A  text  message  containing  the  required  documentation 
appears  in  a  "viewport."  This  documentation  is  available  for  all  SPW  library  blocks  in  three 
forms:  1)  help  on  the  block  usage,  2)  help  on  the  block’s  input  and  output,  and  3)  help  on  the 
block’s  parameters.  The  first  type  of  help  is  obtained  by  pointing  to  the  block  itself  after 
selecting  HELP  on  the  EZ-Menu.  The  second  type  is  obtained  by  pointing  to  any  input  or 
output  port  of  the  block.  The  third  type  is  obtained  by  pushing  into  the  block  and  pointing  to 
any  parameter.  On-line  help  has  also  been  provided  for  all  blocks  designed  and  implemented 
for  use  in  the  two  simulations  presented  in  the  following  sections2. 


2  A  complete  listing  of  the  help  screens  for  each  user  block  used  to  implement  the  simulations  discussed  herein 
is  provided  in  Volume  II  of  this  report.  For  these  SPW  blocks  all  three  types  of  help  screens  described  above  have 
been  combined  into  the  block  usage  help  screen. 
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3.2  The  Conventional  QPSK  Simulation 

The  purpose  of  this  section  is  to  provide  an  overview  of  the  SPW  simulation  of  a 
conventional  QPSK  communications  system.  An  overview  of  the  conventional  QPSK  baseband 
simulation  is  shown  in  Fig.  (3-3).  The  simulation  contains:  a)  a  standard  QPSK  modulator 
which  transmits  a  two-dimensional  (I  and  Q)  signal,  b)  a  channel  which  contains  multiple 
Continuous  Wave  (CW)  jammers,  a  Gaussian  Partial  Band  (PB)  jammer,  and  background 
thermal  (Gaussian)  noise,  and  c)  either  the  standard  linear  receiver  implemented  via  a  two- 
dimensional  correlator  or  one  of  the  memoryless  LO  processor  implementations  followed  by  the 
correlator. 

The  components  of  a  typical  memoryless  LO  processor  block  are  shown  in  Fig.  (3-4). 
The  first  operation  performed  in  the  LO  processor  block  is  rectangular  to  polar  conversion  which 
provides  the  magnitude  and  phase  of  the  received  signal.  Next,  the  corresponding  approximation 
(eg.  traditional  Histogram,  EBH,  FSA,  etc.)  of  the  received  magnitude  PDF  is  constructed. 
Through  the  use  of  the  PDF  approximation,  in  conjunction  with  the  received  magnitude,  the  LO 
nonlinearity  is  computed.  Finally,  the  values  for  the  LO  nonlinear  transform  and  the  phase  of 
the  received  signal  are  converted  back  to  rectangular  form.  Thus,  the  output  of  the  LO 
processor  block  is  (g,,gQ)  =(g(rk)co$dk,g(rl)smdk) ,  where  #(•)  is  the  LO  MNT  estimate. 


Using  (g,  ,gQ),  the  two-dimensional  correlator  generates  the  LO  decision  statistic  of  Eq. 
(1-5)  discussed  in  Section  1.1,  repeated  here  for  convenience: 

N 

{\5(^)cos0t^o^(r,)sin^} 

k»\ 


(3-1) 


Conventional  QPSK  Communications  System 


Figure  (3-3) 


LO  Transform  Block 


Figure  (3-4) 
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3.3  The  QPSK  DSSS  Simulation 


The  simulation  of  a  baseband  QPSK  DSSS  communication  system  with  LO  processing 
has  been  implemented  using  the  SPW  software  package.  An  overview  of  the  entire  system  is 

provided  in  Fig.  (3-5).  It  consists  of:  a)  a  transmitter,  with  output  (X;,XC),  containing  a  QPSK 

source,  whose  output  is  given  by  ( s, ,  sQ ) ,  modulated  by  a  Pseudo-Noise  (PN)  sequence 

generator  producing  a  code  sequence  vector  c ,  b)  a  channel  consisting  of  additive  thermal  (white 

% 

Gaussian)  noise,  multiple  CW  jammers,  and  a  Gaussian  PB  jammer,  and  c)  a  receiver,  with 
input  (rj,~rQ) ,  that  includes  the  LO  MNT  which  outputs  the  signal  vector  the  PN 

sequence  despreading  source  (perfect  synchronization  of  transmitter  and  receiver  is  assumed) 
which  when  multiplied  with  the  output  of  the  LO  MNT  produces  the  ordered  vector  pair  , 
and  a  two-dimensional  correlator. 

LO  QPSK  DSSS  Simulation 


Transmi tter  Channel  Receiver 


Figure  (3-5) 
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3.3.1 


Theoretical  Justification 


To  construct  the  LO  receiver  for  the  QPSK  DSSS  system  it  is  necessary  to  implement 
the  LO  decision  statistic  of  Eq.  (1-5),  repeated  here  for  convenience 


N 

(,=E  {sLg(rt)cosdk+sQ J(rk)sind\ 

*•1  v 


(3-2) 


Referring  to  Fig.  (3-5),  the  k*  sample  of  the  baseband  output  of  the  QPSK  DSSS  transmitter  is 

% 

given  by  (X,^,X0  ) -(cksL,cksQJ ,  where  c=[c,,c2>...,cJt,...,cJJr  is  the  PN  sequence  which  is  fixed 
for  a  given  signal  period,  and  (s^,sQ  )  is  the  k*  sample  of  the  m*  possible  transmitted  signal 
pair.  Since  ck  •  (cks^,cksQ  )  =  (s,  ,SqJ  for  the  DSSS  system,  there  is  a  one-to-one  and  invertible 
mapping  of  s,  to  X,  .  The  process  of  memoryless  LO  detection  of  the  received  signal  can  be 
implemented  by  correlating  the  output  of  the  LO  MNT,  over  the  set  of 

{(X'.,x<?.)}>  m=l,...,Af,  choosing  the  one  corresponding  to  the  largest  decision  statistic,  l„,  and 

mapping  the  signal  vector  (X,  ,\Q )  to  the  possible  information  signal  vector  (*£  ,~sQ ) .  This 
procedure  may  be  stated  more  concisely  as: 

Choose  the  information  signal  (^  ,  sQ  )  corresponding  to  the  signal  vector  (X,  ,XQ )  which 
maximizes: 


/ V 

S(r*)  c°s0* +Xfl  ^(r*)  sin^J 

*-i 


(3-3) 


The  receiver  implementation  of  Eq.  (3-3)  is  termed  the  canonical  form  LO  receiver  since  it  has 
the  same  form  as  Eq.  (3-2).  The  block  diagram  for  the  LO  QPSK  DSSS  receiver 
implementation  given  by  Eq.  (3-3)  is  shown  in  Fig.  (3-6). 
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Canonical  Form  of  the  LO  QPSK  DSSS  System 


Trans m i Channe -  Rece i ver 


Figure  (3-6) 


An  alternate  LO  receiver  implementation  requiring  only  a  slight  modification  to  the 
standard  correlator-based  linear  receiver  can  be  developed  by  noting  that  Eq.  (3-3)  may  be 
written  as 


L  =  E  K  8(rk)  cos  0k  *\j{rk)  sin^} 

N 

58  E  {c*V<r*)c0s^  *ctSaJ(rk)sm9t}  (3-4) 

,V 

= E  {•y^[c^(r*)cos^]+5o.[c^(r*>sin^]} 

Thus,  the  LO  QPSK  DSSS  receiver  can  be  implemented  by  first  multiplying  the  output  of  the 
LO  MNT  by  the  PN  sequence,  and  then  correlating  over  the  set  of  possible  information  signals 

{(T ’*(?>}«  w=l,...,M.  Note  that  this  LO  receiver  may  be  easily  implemented  by  placing  the 

LO  MNT  before  the  despreading  in  a  linear  QPSK  DSSS  receiver.  The  resulting  alternate,  or 
non-canonical  form ,  LO  QPSK  DSSS  receiver  implementation  can  be  summarized  in  the 
following  algorithm: 

Choose  the  information  signal  (s,,sQ)  which  maximizes 
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(3-5) 


*«.  ck8(ri)  cos®*  ck  g(rk)  sind*} 

Due  to  its  modular  nature,  the  LO  receiver  algorithm  given  by  Eq.  (3-5)  was  implemented  in 
the  SPW  QPSK  DSSS  simulation  shown  in  Fig.  (3-5). 

The  non-canonical  form  LO  receiver  algorithm  of  Eq.  (3-5)  has  two  major  advantages 
over  the  canonical  form  algorithm  of  Eq.  (3-3).  First,  it  has  a  more  modular  structure  than  the 
canonical  algorithm  and  can  be  easily  implemented  by  inserting  the  LO  MNT  into  a  standard 
correlator  based  QPSK  DSSS  linear  receiver.  Second,  it  is  more  flexible  than  the  canonical 
algorithm  in  the  sense  that  the  PN  spreading  sequence  can  be  modified  without  requiring  new 
reference  signals  to  be  stored  for  use  by  the  correlator.  However,  the  canonical  algorithm  of 
Eq.  (3-3)  has  the  advantage  of  being  a  potentially  faster  implementation  than  the  non-canonical 
algorithm,  since  only  one  multiplication  is  required  (i.e.  multiplying  the  output  of  the  LO  MNT 
by  the  stored  reference  signals),  whereas  the  non-canonical  algorithm  requires  two 
multiplications  (i.e.  multiplication  of  the  output  of  the  LO  MNT  by  the  PN  sequence,  and  then 
multiplying  the  result  by  the  stored  reference  signals).  System  requirements,  e.g.  speed  versus 
flexibility,  will  dictate  which  LO  receiver  implementation  is  more  desirable. 
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4. 


SIMULATION  RESULTS  FOR  THE  MEMORYLESS  LOCALLY  OPTIMAL 
PROCESSING  ALGORITHMS 


This  section  presents  the  simulation  results  for  the  conventional  QPSK  and  DSSS  QPSK 
systems  with  memoryless  LO  processing  discussed  in  Sections  3.2  and  3.3.  The  following 
memoryless  LO  algorithms  were  examined  in  each  system:  Histogram,  EBH,  and  FSA.  In 
addition,  the  second  order  M-Interval  Polynomial  Approximation  (MIPA)  algorithm  [Illi91, 
Illi93]  was  used  in  conjunction  with  the  DSSS  QPSK  system.  The  probability  of  bit  error,  Pbt 
performance  of  each  system  and  LO  algorithm  combination  was  computed  for  a  multitude  of 
jamming  scenarios.  When  useful,  the  Pb  for  a  standard  linear  system  was  also  calculated  for  the 
purpose  of  comparison.  By  examining  the  operation  of  the  LO  algorithms  in  many  different 
environments,  it  is  possible  to  determine  basic  characteristics  which  govern  performance,  and 
ultimately  which  system/LO  algorithm  configuration  is  optimal  for  a  given  jamming  scenario. 
This  section  is  organized  in  the  following  manner.  First,  the  simulation  results  for  each  LO 
algorithm  in  the  conventional  QPSK  system  are  addressed.  Next,  the  results  are  presented  for 
the  QPSK  DSSS  system  when  subjected  to  three  different  types  of  jammers,  namely  a  single  CW 
jammer,  two  CW  jammers,  and  a  PB  jammer.  Finally,  the  performance  of  both  systems  are 
compared  to  determine  the  effects  of  the  SS  system’s  processing  gain  on  the  LO  processor. 


4.1  Results  for  the  Conventional  QPSK  System  -  The  Single  CW  Jammer  Scenario 

The  results  for  the  conventional  QPSK  system  when  subjected  to  a  single  CW  jammer 
are  presented  in  this  section.  The  user-definable  parameters  for  each  LO  algorithm  are  shown 
in  Table  4-1.  As  may  be  seen,  the  parameters  B  (number  of  bins),  R,  (Symbol  Rate)1,  N, 
(Samples  per  Symbol),  and  N  (Samples  per  Correlation)  are  common  to  all  the  memoryless  LO 
algorithms. 


1  AH  other  simulation  parameters  of  interest  are  normalized  to  R,.  Thus,  the  value  chosen  for  R,  is  arbitrary. 
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^aramete^^o?th?T3^1gont5nM 


LO  Algorithm 

1) 

Histogram 

Symbol  Rate 

2) 

EBH 

B  - 

Number  of  Bins 

3) 

MIPA  (2nd  Order) 

K’ 

Samples  per  Symbol 

N  - 

Samples  per  Correlation  (equivalent 
to  samples  per  PDF  approximation) 

4) 

FSA 

Symbol  Rate 

P  - 

Order  of  Approximation 

B  - 

Number  of  Bins 

Samples  per  Symbol 

N  - 

Samples  per  Correlation  (equivalent 

to  samples  per  PDF  approximation) 


Table  4-1 

Before  it  is  possible  to  thoroughly  test  each  LO  algorithm  for  a  variety  of  Jammer-to- 
Signal  power  ratio  (7/5),  CW  jammer  frequency  relative  to  the  symbol  rate  if/Rs),  and  bit  energy 
to  signal  thermal  (Gaussian)  background  noise  power  ratio  (E,/N0)  values2,  it  was  necessary  to 
determine  values  for  B,  N„  and  N.  In  addition,  P  (Order  of  Approximation)  must  be  chosen  for 
the  FSA  LO  algorithm.  To  determine  a  value  of  N,  that  would  be  used  in  all  subsequent 
analyses,  the  Pb  for  both  the  Histogram  and  EBH  methods  were  computed  using  the  system 
settings  shown  in  Table  4-2,  and  the  results  are  presented  in  Fig.  (4-1),  with  the  region 
15  <  Ns  £  50  emphasized  in  Fig.  (4-2).  Note  that  the  rows  are  reversed  in  Fig.  (4-2)  for  clarity. 
As  may  be  seen  from  the  figures,  Pb  for  the  conventional  QPSK  system  with  LO  processing 
decreases  as  Ns  increases.  Using  these  results,  and  accounting  for  processing  limitations,  N,  = 
20  was  used  for  all  subsequent  simulations. 


2  In  this  report,  the  values  for  J/S  and  EJN0  are  specified  per  channel. 
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Simulation  Parameters 


Parameter 

Value 

N, 

5  to  50 

R, 

0.1 

B 

64 

N 

50,400 

EM 

10  dB 

N0 

0  dB 

J/S 

30  dB 

m 

0.496 

Symbols  per 

7,056 

Pb  Calculation 

Table  4-2 


To  choose  values  for  B  and  N,  the  Pb  was  computed  for 
the  Histogram  and  EBH  methods  for  a  range  of  B  and  N  values 
when  the  system  was  subjected  to  a  single  CW  jammer  with  J/S 
=  30  dB.  The  values  for  the  system  parameters  used  in  this 
analysis  are  given  in  Table  4-3,  and  the  results  for  the 
Histogram  and  EBH  methods  are  displayed  in  Figs.  (4-3)  and 
(4-4),  respectively.  The  Pb  plots  indicate  that  there  is  a  critical 
relationship  between  B  and  N  for  both  methods.  If  the  chosen 
value  of  B  is  too  large  for  a  given  value  of  N  the  Pb  suffers. 
Conversely,  if  the  chosen  value  of  N  is  too  small  for  a  given 
value  of  B,  the  Pb  also  increases. 


Conventional  QPSK  System 
Histogram  and  EBH  Implementations 


Conventional  QPSK  System 
Histogram  and  EBH  Implementations 


Figure  (4-1) 


Figure  (4-2) 
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The  main  reason  for  these  relationships  is  that  in  both  the  Histogram  and  EBH  methods, 
the  total  number  of  samples,  N,  is  separated  into  B  bins.  As  B  increases  for  a  given  N,  the 
average  number  of  samples  in  each  bin  decreases.  If  this  average  number  becomes  too  small, 
the  PDF  approximation  will  not  be  a  good  estimate  of  the  actual  PDF,  and  the  resulting  MNT 
will  be  inaccurate.  For  all  subsequent  analyses ,  the  value  of  N  was  set  to  N  —  50,000  and  the 
effect  of  varying  B  on  the  Pb  was  examined. 

With  the  basic  analysis  of  the  LO  processor  parameters  completed,  it  was  possible  to 
examine  its  performance  relative  to  a  variety  of  channel  parameters,  including:  J/S ,  f/R„  and 
EJN0.  In  the  first  simulation,  the  relationships  between  these  three  parameters  and  B  were 
examined  for  the  Histogram  and  EBH  methods.  For  the  system  parameters  in  Table  4-4,  the 
following  analyses  are  plotted:  1)  for  E,/N0  =  0  dB,  Figs.  (4-5)  and  (4-6)  show  the  Pb  relative 
to  J/S  and  f/R ,  3  for  5  =  8,  Figs.  (4-7)  and  (4-8)  show  the  Pb  relative  to  J/S  and  f/Rs  for  B  = 
32,  and  Figs.  (4-9)  and  (4-10)  show  Pb  relative  to  J/S  and  f/Rs  for  B  =  128;  2)  for  E,/N0  =  10 
dB,  Figs.  (4-11)  and  (4-12)  show  the  Pb  relative  to  J/S  and  f/R,  for  5  =  8,  Figs.  (4-13)  and 
(4-14)  show  the  Pb  relative  to  J/S  and  f/R,  for  B  =  32,  and  Figs.  (4-15)  and  (4-16)  show  Pb 
relative  to  J/S  and  f/R ,  for  B  =  128.  For  a  baseline  measurement,  a  plot  of  Pb  for  a 
conventional  linear  QPSK  system  with  Eb/N0  =  10  dB  is  provided  in  Fig.  (4-17). 

The  first  observation  that  can  be  made  from  these  figures  is  that  there  is  a  clearly  defined 
region  in  which  the  LO  processor  improves  performance.  In  all  the  figures,  this  region  begins 
around  J/S  >  20  dB  and  f/R,  >  0.2.  Note  that  the  relationship  between  Pb  and  f/R,  implies 
that  the  performance  of  the  LO  processor  improves  as  the  CW  jammer  moves  farther  away  from 
the  center  of  the  main  lobe  bandwidth  of  the  information  signal,  and  only  a  slight  performance 
improvement  is  seen  as  the  CW  jammer  moves  towards  the  center  of  the  main  lobe.  The  next 


3  The  range  of  f/R,  was  chosen  to  be  0.00496  to  0.496,  rather  than  0.005  to  0.5,  so  that  samples  of  the  CW 
jammer  taken  in  successive  periods  would  not  have  the  same  value.  This  produces  a  histogram  that  has  a  relatively 
'smooth*  shape,  rather  than  one  that  has  a  few  bins  of  high  probability,  and  a  majority  with  zero  probability. 


4-5 


observation  is  that  increasing  B  only  decreases  Pb  up  to  a  point,  after  which  the  Pb  increases. 
Furthermore,  the  effects  of  varying  B  seem  to  be  more  pronounced  at  low  E,JN0  than  at  high 
EJNq.  This  can  be  seen  more  clearly  in  Fig.  (4-18)  to  (4-21)  where  Pb  is  shown  for  B  and  f/R} 
varying. 

Another  interesting  phenomenon  that  can  be  observed  is  that  while  J/S  must  be  large  for 
good  performance,  the  LO  processor  is  able  to  reduce  the  Pb  more  substantially  in  a  high  E^/Nq 
environment  (e.g.  10  dB)  than  in  a  low  EJN0  environment  (e.g.  0  dB).  This  indicates  that  for 
the  best  LO  processor  performance,  the  signal  must  be  much  weaker  than  the  CW  jammer,  but 
at  the  same  time  larger  than  the  background  thermal  (white  Gaussian)  noise. 


Simulation  Parameters 


Parameter 

J/S 

m* 

EM 

B 


Value 

0  to  SO  dB 
0.00496  to 
0.496 
0  to  10  dB 
8  to  128 


R,  0.05 

Ns  20 

N  50,000 

N0  0  dB 

Symbols  per  5,000 


Pb  Calculation 


Table  4-4 


Finally,  in  comparing  the  results  for  the 
Histogram  method  to  the  EBH  method  for  the 
parameters  in  Table  4-4,  it  may  be  observed  that 
the  Histogram  method  has  a  slightly  lower  Pb  than 
the  EBH  method  for  most  cases  in  the  region  of 
interest.  It  must  be  stressed,  however,  that  the 
value  of  Ns  was  N,  =  20  for  all  simulation  runs  in 
Figs.  (4-5)  to  (4-16).  But  from  Fig.  (4-2)  it  may 
be  seen  that  while  the  Histogram  method  performs 
better  than  the  EBH  method  for  low  values  of  Ns, 
the  EBH  method  produces  a  lower  Pb  for  higher 
values  of  N}.  Thus,  by  increasing  Ns  to  a  value 
larger  than  20,  it  may  be  possible  to  improve  the 
performance  of  the  EBH  method  such  that  it  will 
have  a  lower  Pb  than  the  Histogram  method  for 


most  cases  in  the  region  of  interest. 
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Probability  of  Bit  Error 


Conventional  QPSK  System  Conventional  QPSK  System 

Histogram  Method  -  Eb/N0=0  dB,  B=8  EBH  Method  •  Eb/N0»0  dB,  B*8 


Figure  (4-5)  Figure  (4-6) 


Conventional  QPSK  System  Conventional  QPSK  System 

Histogram  Method  -  Eb/N0=0  dB,  B=32  EBH  Method  -  Eb/N0*0  dB,  B=*32 


J/S(dB)  1/S(d8) 


Figure  (4-7)  Figure  (4-8) 
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Probability  o{  Bit  Error  Probability  ol  Bit  Error 


Conventional  QPSK  System 
Histogram  Method  -  Eb/N0=0  dB,  B=128 


W 


Conventional  QPSK  System 
EBH  Method  •  Eb/N0-0  dB,  B-128 
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Figure  (4-9) 

Figure  (4-10) 

Conventional  QPSK  System 

Histogram  Method  -  Eb/N0»10  dB,  B=8 

Conventional  QPSK  System 
EBH  Method  -  Eb/N0=10  dB,  B*8 
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Figure  (4-11) 


Figure  (4-12) 


Probability  of  Bit  Ettot 


Conventional  QPSK  System 
Histogram  Method  -  Eb/N0*10  dB,  B*32 


Figure  (4-13) 


Conventional  QPSK  System 
EBH  Method  -  Eb/N0=10  dB,  B-32 


Figure  (4-14) 


Conventional  QPSK  System 
Histogram  Method  -  Eb/N0=l0  dB,  B— 128 


Conventional  QPSK  System 
EBH  Method  -  Eb/N0»10  dB,  B»128 


J/S  (dB) 


J/S  (dB) 


Figure  (4-15) 


Figure  (4-16) 


4-9 


Conventional  QPSK  System 
Linear  Receiver  -  Eb/N0=  10  dB 


Conventional  QPSK  System 
Histogram  Method  -  Eb/N0=»10dB,  J/S 
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Figure  (4-17) 

Figure  (4-18) 

Conventional  QPSK  System 

EBH  Method  -  Eb/NO-lOdB,  J/S»0dB 

Conventional  QPSK  System 
Histogram  Method  -  Eb/N0»10dB,  J/S*50dB 
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Figure  (4-19) 


Figure  (4-20) 


Probability  of  Bit  error 


The  Pb  performance  for  the  histogram 
implementation  of  the  FSA  LO  algorithm  in 
conjunction  with  the  conventional  QPSK  system 
was  also  investigated.  Before  examining  the  Pb 
performance  with  respect  to  J/S,f/R„  and  Eb/N0, 
it  was  necessary  to  determine  "good"  values  for 
P  and  B.  Figure  (4-22)  shows  the  Pb 
performance  of  the  FSA  method  relative  to  P  and 
B  for  a  single  CW  jammer  having  J/S  =  30  dB 
and  f/R,  —  0.496.  The  system  parameters  are 
displayed  in  Table  4-5.  As  may  be  seen  from 
Fig.  (4-22),  there  is  a  definite  region  where  the 
FSA  order,  P,  exceeds  the  histogram  resolution, 
or  B.  Using  these  results,  and  after  further 
investigation,  the  values  P  =  8  and  B  =  64  were  chosen  for  the  subsequent  analysis. 

Using  the  knowledge  gained  in  the  preceding  analysis,  the  Pb  for  the  histogram 
implementation  of  the  FSA  method  was  computed  relative  to  7/S,  f/R„  and  EJNo.  The  results 
are  displayed  in  Figs.  (4-23)  to  (4-25)  for  the  simulation  parameters  shown  in  Table  4-6.  The 
first  major  observation  is  that  the  FSA  method  is  governed  by  the  same  region  of  interest  as  the 
Histogram  and  EBH  methods,  namely  J/S  >  20  dB  and  f/Rs  >  0.2,  and  the  Pb  decreases  as 
f/Rs  increases  as  before.  Since  these  two  observations  are  the  same  as  those  seen  in  the 
Histogram  and  EBH  methods,  it  may  be  assumed  they  are  characteristic  of  the  basic  LO 
processor  algorithm  and  not  the  result  of  a  specific  implementation.  One  final  observation  of 
interest  is  that  Pb  decreases  as  EJN0  increases.  Thus,  while  it  is  critical  that  J/S  be  large  enough 
for  the  small-signal  assumption  to  be  valid,  it  is  also  important  that  the  information  signal  be 
larger  than  the  background  thermal  (white  Gaussian)  noise. 


Conventional  QPSK  System 
EBH  Method  -  Eb/N0=l0dB,  J/S=50dB 


•  t*  si  m  ixa 


Number  of  Bias 

Figure  (4-21) 
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™l^""^l^SulatIo^Sra5etere™,,™ 
Parameter  Value 


p 

2  to  20 

B 

8  to  128 

f/R, 

0.00496  to 

0.496 

K 

0.05 

N, 

20 

N 

50,000 

EM 

10  dB 

J/S 

30  dB 

No 

0  dB 

Symbols  per 

Pb  Calculation 

5,000 

Conventional  QPSK  System 
FSA  Method  -  J/S*30  dB,  fj/Rs  *0.496 


Table  4-5 


Figure  (4-22) 


In  comparing  all  three  methods,  the  Histogram,  EBH,  and  the  histogram  implementation 
of  the  FSA,  it  may  be  observed  that  all  three  have  approximately  the  same  Pb  performance. 
This  is  most  likely  a  result  of  each  method  being  based  on  a  histogram  approximation  of  the 
PDF.  Because  the  parametric  implementation  of  the  FSA  is  not  based  on  an  initial  histogram 
approximation,  it  is  believed  that  this  method  will  have  a  lower  Pb  in  the  region  of  interest  than 
the  other  LO  processor  algorithms.  This  claim  will  be  examined  further  in  future  research. 
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Probability  of  BU  Error 


4.2  Results  for  the  QPSK  DSSS  System 


The  performance  of  the  QPSK  DSSS  system  was  examined  for  the  following  overall 
jamming  scenarios:  1)  a  single  CW  jammer,  2)  two  CW  jammers,  and  3)  a  single  PB  jammer. 

4.2.1  The  Single  CW  Jammer  Scenario 

The  results  for  the  QPSK  DSSS  system  when  subjected  to  a  single  CW  jammer  are 
presented  in  this  section.  Except  for  the  addition  of  a  new  LO  processor  parameter,  the 
Processing  Gain  (PG),  the  channel  parameters  (7/5,  f/Rs,  and  EJN/)  and  the  LO  processor 
parameters  (Rlt  N„  N,  B,  and  P)  of  interest  are  the  same  as  for  the  conventional  QPSK  system. 
In  this  report,  PG  is  defined  as  the  number  of  chips  per  information  symbol,  or  in  terms  of  the 
other  LO  parameters: 

N 

PG  =  -l  (4-D 

Nc 

where  Ne  is  the  number  of  samples  per  chip. 

The  Pb  performance  of  the  QPSK  DSSS  system  with  LO  processing  relative  to  J/S,f/Rs, 
and  PG  was  examined  using  the  system  parameters  given  in  Table  4-7.  The  values  for  Nt,  N, 
B,  and  P  were  chosen  using  the  knowledge  gained  through  the  analysis  of  the  conventional 
QPSK  system.  The  organization  of  the  Pb  plots  is  as  follows.  For  the  Histogram  method.  Figs. 
(4-26)  and  (4-27)  show  the  Pb  with  respect  to  7/5  and  f/R, 4  and  Figs.  (4-28)  and  (4-29)  show 
the  Pb  with  respect  to  PG  and  f/Rs.  Similarly,  Figs.  (4-30)  to  (4-33),  (4-34)  to  (4-37),  and 

4  Note  that  in  these  analyses  Pb  is  plotted  versus  f/R,  for  f/R,  in  the  range  0.201376  to  0.496,  not  0.00496  to 
0.496  as  in  the  results  for  the  conventional  QPSK  system.  Values  of  f/R,  less  than  0.2  are  outside  the  region  of 
interest  because  the  Pb  is  large  there. 
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(4-38)  to  (4-41)  depict  Pb  with  respect  to  J/S,  PG,  and  f/R,  for  the  EBH,  FSA,  and  second-order 
MIPA  (MIPA2)  methods,  respectively.  Comparable  results  for  a  standard  linear  QPSK  DSSS 
system  are  provided  in  Figs.  (4-42)  to  (4-45).  From  a  cursory  examination  of  the  Pb  plots,  three 
observations  can  be  made:  1)  for  f/R,  >  0.2  and  a  given  PG,  the  Pb  for  almost  all  LO 
processor  methods  decreases  as  J/S  increases;  2)  for  J/S  >  10  dB  and  a  given  PG,  the  Pb 
decreases  as  f/R,  increases;  and  3)  for  a  given  J/S  and  f/R„  the  Pb  for  almost  all  LO  processor 
methods  decreases  as  PG  increases.  These  results  are  consistent  with  those  for  the  conventional 
QPSK  system. 


SimuSationTarametere 


Parameter 

Value 

J/S 

0  to  50  dB 

f/K 

0.201376 
to  0.496 

PG 

5  to  20 

R, 

0.1 

K 

2 

K 

PG  *  Nc 

N 

50,400 

B  (histogram  and  EBH) 

32 

(FSA) 

64 

(MIPA2) 

8 

P  (FSA) 

8 

EM 

10  dB 

No 

0  dB 

Symbols  per 

Pb  Calculation 

5,040 

Table  4-7 


In  comparing  across  LO  processor 
methods,  the  following  observations  can  be  made. 
In  the  region  of  interest,  J/S  >  20  dB  and  f/R,  > 
0.2,  all  the  LO  processor  methods  have  roughly 
the  same  Pb  performance  characteristics  with  the 
following  exceptions:  1)  the  FSA  and  MIPA2 
methods  perform  better  than  the  Histogram  and 
EBH  methods  for  low  J/S  ( J/S  <  10  dB),  and  2) 
the  Histogram,  FSA,  and  MIPA2  methods  perform 
better  than  the  EBH  method  when  PG  and  J/S  are 
both  low  (PG  ^  10,  J/S  <  10  dB).  The  second 
exception  is  consistent  with  the  results  observed 
for  the  conventional  QPSK  system,  namely  that  for 
low  N,  the  Histogram  method  has  a  smaller  Pb 
than  the  EBH  method,  but  the  converse  is  true  for 
high  N,.  Future  research  is  required  to  determine 
the  reason(s)  for  the  first  exception. 
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Probability  of  Bit  iirror  Probability  of  Bit  Error 


Probability  of  Bit  Ertor  Probability  of  Bit  Error 


4.2.2  The  Two  CW  Jammers  Scenario 

This  section  presents  the  results  for  the  QPSK  DSSS  system  using  the  Histogram  and 
EBH  implementations  of  the  LO  nonlinear  processor  when  subjected  to  two  CW  jammers.  In 
this  analysis,  the  first  CW  jammer  was  held  fixed  at  J/S  (1*  CW)  =  30  dB  and  f/R,  (1*  CW) 
=  0.496.  The  second  CW  jammer  was  then  varied  over  the  range  J/S  (2nd  CW)  =  -20  to  50 
dB  and  f/Rs  (2nd  CW)  =  0.00406  to  0.406.  The  remaining  system  parameters  are  given  in  Table 
4-8.  The  Pb  performance  curves  are  plotted  in  Figs.  (4-46)  and  (4-47)  for  the  Histogram  and 
EBH  methods,  respectively.  As  can  be  seen  from  the  figures,  the  Pb  is  small  when  J/S  (2nd  CW) 
is  low,  and  increases  as  J/S  (2nd  CW)  increases.  Note  however,  as  J/S  (2nd  CW)  grows  larger 
than  J/S  (1“  CW),  i.e.  as  the  second  jammer  becomes  the  dominant  interferer,  the  Pb  remains 
at  a  high  value.  This  indicates  that  the  absolute  values  of  the  J/S  for  both  CWs,  as  well  as  the 
difference  in  their  magnitudes,  determines  the  Pb  performance. 

Another  observation  that  can  be  made  in  this  scenario  is  that  as  f/Rs  (2nd  CW)  decreases, 
the  Pb  actually  decreases.  This  would  seem  to  contradict  the  results  of  the  previous  sections. 
However,  a  possible  reason  for  this  seemingly  anomalous  case  is  that  as  f/R,  (2nd  CW) 
decreases,  the  second  CW  jammer’s  frequency  spectrum  moves  farther  away  from  that  of  the 
first  CW  jammer.  Under  this  hypothesis  it  appears  that  the  increased  frequency  spacing  between 
interferes  is  not  as  detrimental  to  the  system  as  are  closely  spaced  interferes. 
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Parameter 
J/S  (2nd  CW) 
f/R,  (2nd  CW) 


value 


-20  to 
50  dB 
0.00406  to 
0.406 

0.1 

2 

Ne  *  PG 
20 

50,400 

32 


E/N0  10  dB 

N0  0  dB 

J/S  (1-  CW)  30  dB 

f/R,  (1*  CW)  0.496 

Symbols  per  5,040 


Pb  Calculation 


Table  4-8 


Histogram  QPSK  DSSS  System 
Probability  of  Bit  Error 


J/S  (dB)  for  Second  CW  Jammer 

Figure  (4-46) 


EBH  QPSK  DSSS  System 
Probability  of  Bit  Error 


-20  -10  0  10  20  30  40  SO 

J/S  (dB)  for  Second  CW  Jammer 


Figure  (4-47) 


4.2.3  The  PB  Jammer  Scenario 

This  section  discusses  the  Pb  performance  of  the  QPSK  DSSS  system  when  the  channel 
interference  is  a  single  PB  jammer,  implemented  by  passing  white  Gaussian  noise  through  a 
lowpass  filter,  and  white  Gaussian  background  noise.  The  results  shown  are  for  the  LO 
processor  implemented  via  the  Histogram  method.  Two  parameters  govern  the  PB  jammer:  J/S 
and  /.JR,,  where  fc  is  the  lowpass  bandwidth  (cutoff  frequency)  of  the  PB  jammer  (recall  that  the 
simulation  is  performed  at  baseband).  The  simulation  results  for  f^R,  <  1,  with  the  simulation 
parameters  in  Table  4-9,  are  shown  in  Figs.  (4-48)  to  (4-50).  Figure  (4-48)  illustrates  the 
relationship  of  Pb  to  J/S  and  fJRs  5  for  PG  =  20  (similar  results  were  observed  for  PG  *  5  to 
15  but  are  not  shown  here).  The  plots  of  Pb  with  respect  to  PG  and  f/R,  for  J/S  =  -10  and  30 
dB  are  shown  in  Figs.  (4-49)  and  (4-50),  respectively.  As  may  be  seen,  Pb  is  small  for  low 
values  of  J/S  ( J/S  <0  dB)  and  increases  as  J/S  increases.  In  addition,  it  may  be  observed  that 
Pb  decreases  as  PG  increases  for  high  J/S,  but  for  low  J/S  Pb  actually  increases  as  PG  increases. 
Future  analysis  is  required  to  determine  the  reason(s)  behind  this  trend.  Finally,  for  a  given  PG 
and  J/S,  the  Pb  is  seen  to  increase  as  fJRs  increases.  This  indicates  that  the  LO  processor 
performs  well  for  narrowband  interference  and  worse  for  wideband  interference. 

Figures  (4-51)  to  (4-53)  depict  the  Pb  plots  for  fJRt  >  1  with  the  system  parameters 
given  in  Table  4-10.  These  results  confirm  that  the  performance  of  the  LO  processor  degrades 
as  the  power  spectrum  of  the  PB  jammer  becomes  more  wideband. 


3  Unlike  the  CW  jammer,  the  values  for  fJR,  for  the  PB  jammer  may  be  chosen  to  be  "nice"  values,  such  as 
0.2,  0.8,  1.2,  etc.  Since  the  PB  jammer  is  implemented  by  passing  white  Gaussian  noise  through  a  lowpass  filter, 
the  probability  of  successive  samples  being  identical  is  small.  Thus,  any  reasonable  value  of  fJR,  will  produce  a 
histogram  that  is  relatively  'smooth.  “ 
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rameters 


Parameter 

Value 

J/S 

-10  to 

30  dB 

fJK 

0.2to0.8 

R, 

0.1 

Ne 

2 

Ns 

Ne  *  PG 

PG 

20 

N 

50,400 

B 

32 

EM 

10  dB 

N0 

0  dB 

Symbols  per 

Pb  Calculation 

5,040 

Table  4-9 


Histogram  QPSK  DSSS  System 
PB  Jammer  -  J/S=-lQ  dB 
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Pmanunf  Gtia 

Figure  (4-49) 


Histogram  QPSK  DSSS  System 
PB  Jammer  -  PG  »  20 


t#  •  i#  » 

J/S(dB) 


Figure  (4-48) 


Histogram  QPSK  DSSS  System 
PB  Jammer  -  J/S»30  dB 
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Si 
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Processing  Gaia 


Figure  (4-50) 


4.3  Comparison  of  the  Conventional  QPSK  and  QPSK  DSSS  Systems 

To  compare  the  effects  of  spread  spectrum  on  the  LO  processor,  the  Pb  performance  of 
the  conventional  QPSK  and  QPSK  DSSS  systems  using  the  Histogram,  EBH,  and  FSA  LO 
implementation  methods  were  compared  for  the  single  CW  jammer  scenario.  As  in  previous 
simulations,  N,  for  the  conventional  QPSK  system  was  fixed  at  N,  —  20,  and  PG  =  10  and  Nc 
-  2  were  used  for  the  QPSK  DSSS  system.  Figs.  (4-54)  through  (4-59)  illustrate  the  Pb  versus 
J/S  results  with  f/Rs  -  0.397792  and  0.496  for  each  of  the  three  LO  processor  methods.  The 
corresponding  system  parameters  are  provided  in  Table  4-11. 


Simulation  Parameters 


As  may  be  seen  in  the  figures,  the  Pb  for 


Parameter 

Yfltw 

J/S 

0  to  50  dB 

m 

0.397792 

and  0.496 

R, 

0.1 

N,  (conv.  QPSK) 

20 

Ne  (QPSK  DSSS) 

2 

PG  (QPSK  DSSS) 

10 

N 

50,400 

B  (histogram  and  EBH) 

32 

(FSA) 

64 

P  (FSA) 

8 

EM 

10  dB 

N0 

0  dB 

Symbols  per 

5,040 

Pb  Calculation 

Table  4-11 

the  QPSK  DSSS  system  is  smaller  than  that  of  the 
conventional  QPSK  system  for  low  values  of  J/S. 
However,  for  larger  values  of  J/S  the  Pb  for  the 
two  systems  is  approximately  the  same.  Thus,  for 
high  values  of  J/S ,  the  LO  processor  dominates  the 
improvement  in  performance.  In  addition,  spread 
spectrum  can  be  used  to  improve  performance,  and 
provide  robustness,  at  lower  values  of  J/S. 
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Probability  of  Bit  Error 


Comparison  of  Both  Systems 
Histogram  Method  -  fj/Rs  *0.397792 


•  it  n  m  <•  n 


J/S  (dB) 


Comparison  of  Both  Systems 
EBH  Method  -  fj/Rs =0.397792 


A. 


Figure  (4-54) 


Figure  (4-55) 


Comparison  of  Both  Systems 
FSA  Method  -  fj/Rs  =0.397792 
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Comparison  of  Both  Systems 
Histogram  Method  -  fj/Rs =0.496 
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Figure  (4-56) 


Figure  (4-57) 
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Comparison  of  Both  Systems  Comparison  of  Both  Systems 

EBH  Method  -  fj/Rs =0.496  FSA  Method  -  fj/Rs  *0.496 


J/S(dB)  J/S(d8) 


Figure  (4-58) 


Figure  (4-59) 


4.4  Conclusions 


Using  the  previous  analyses,  a  number  of  conclusions  can  be  drawn  concerning  the 
optimum  use  of  the  memoryless  LO  processor.  When  the  channel  interference  consists  of  a 
single  CW  jammer: 

•  the  memoryless  LO  processor  is  only  effective  in  high  J/S  and  high  f/R, 
environments 

•  the  performance  of  the  LO  processor  improves  as  the  EJNq  increases 

•  large  values  of  N \  should  be  used 

•  regarding  the  LO  processor  implementation  methods,  for  a  given  value  of  N,  the 
value  for  B  should  be  chosen  judiciously  such  that  the  average  number  of  samples 
in  each  histogram  bin  remains  large  enough  so  that  an  accurate  approximation  of 
the  PDF  is  obtained 

In  the  case  of  two  CW  jammers,  the  memoryless  LO  processor  should  be  used  only  if  one 
of  the  jammers  has  a  relatively  low  J/S.  In  addition,  the  Pb  performance  improves  as  the 
spectral  separation  of  the  CW  jammers  increases. 

If  the  channel  contains  a  single  PB  jammer,  the  memoryless  LO  processor  should  only 
be  used  in  moderate  J/S  environments.  Also,  the  PG  of  the  DSSS  system  must  be  modified  as 
the  J/S  changes. 

Finally,  it  should  be  noted  that  all  of  the  results  in  this  section  are  for  the  memoryless 
LO  processor  only.  Future  research  should  examine  the  use  of  LO  processors  with  memory, 
and  their  ability  to  improve  Pb  performance,  particularly  in  regions  where  the  memoryless  LO 
processor  is  insufficient. 
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5. 


LOCALLY  OPTIMAL  PROCESSING  WITH  MEMORY 


The  LO  processors  discussed  in  the  previous  sections  assumed  that  adjacent  samples  of 
the  total  channel  interference  were  Independent  and  Identically  Distributed  (IID).  Using  this 
assumption,  it  is  possible  to  obtain  the  joint  PDF  of  the  interference  samples  as  the  product  of 
the  individual  marginal  PDFs.  However,  the  resulting  memoryless  LO  processor  does  not  take 
advantage  of  the  wealth  of  information  contained  in  the  interrelationships  between  adjacent 
samples.  Thus,  a  more  robust  LO  processor,  of  which  the  memoryless  processor  is  a  subset, 
can  be  derived  by  removing  the  IID  assumption.  This  results  in  a  processor  with  memory.  Of 
particular  interest  during  the  contract  effort  was  the  derivation  of  the  LO  processor  with  memory 
presented  in  the  Robust  Digital  Adaptive  Transceiver  (RDAT)  report  by  Charles  Stark  Draper 
Laboratories,  Inc.  [Drap89],  The  following  sections  present  the  RDAT  theory  and  methodology, 
focusing  on  the  development  of  the  LO  detector  for  one-dimensional  (ID)  and  two-dimensional 
(2D)  signalling  environments.  In  addition,  a  particular  interference  scenario  was  simulated  to 
provide  further  insight  into,  and  a  means  of  verifying,  the  research  presented  in  the  RDAT 
report. 


5.1  Derivation  of  the  LO  Processor  with  Memory  for  ID  Signalling1 

The  LO  processor  with  memory  for  the  case  of  ID  signalling  is  derived  in  [Drap89]  as 
the  solution  to  the  following  radar  problem:  After  translation  to  baseband,  consider  the 
following  two  hypotheses: 


1  The  notation  for  all  of  Section  5  follows  the  symbology  introduced  in  [Drap89J.  When  necessary  to  clarify 
the  meaning  of  an  equation  or  expression,  notation  consistent  with  this  report  is  also  included. 
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Hr  VM  (j.,) 

W0-  *v 

where  "^■[z,,...,zjr  the  vector  °f  received  signal  samples,  is  the  vector  of 

transmitted  signal  samples  and  may  be  written  as  =r[c1,...,cN]r  with  the  set  of  {cj  s 

having  amplitudes  equal  to  i  1  and  where  s  is  a  constant,  is  the  vector  of  total 

channel  noise  (thermal  noise  and  jammer)  samples,  and  N  is  the  number  of  samples  per  symbol 
period.  The  N-point  joint  PDF  of  the  channel  noise  is  given  by  /.(**)•  The  following 

assumptions  are  made  in  the  derivation: 

•  cv  and  (and  thus  %  and  n^)  are  statistically  independent 

•  the  signal  of  interest  is  much  weaker  than  the  total  channel  noise. 

The  LO  processor  is  derived  by  considering  the  standard  likelihood  ratio  test  [VanT68]: 

W, 

, .  /,i«,(whi>  >  (5-2) 

no 


or  equivalently 


,n(o"n  m 
•  ■J 


since  the  natural  logarithm  is  a  monotonic  function.  Using  the  assumption  that  the  signal  and 
interference  are  statistically  independent,  /  becomes 
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(5-4) 


~7m~ 

Taking  the  natural  logarithm  of  /  and  using  the  first  order  N-dimensional  Taylor  series  expansion 
gives  the  expression 


ln(/)  =  In <fn(  ^))  -  In (£(  4)) 


k’WJLV) 


-  ■fcW.M  * 


j 

where  ~-~h( zj  =  JLh(zN),...,-£-h(zN)  is  the  gradient  of  /t( zj .  The  expression  of  Hq.  (5-5) 
uZf,  oz,  dzs 


is  a  valid  approximation  due  to  the  weak  signal  assumption.  Substituting  Eq.  (5-5)  into  Eq. 
(5-3)  results  in  the  following  form  for  the  LO  detector  with  memory: 


or  equivalently 


=sG^>ln(Ti) 


1 

r  „  7*T~T  >  In(-n) 
L^Gscs  <— ^=7 
^  s 


where  GN-[gx,...,gN]T  with  the  nonlinearity  at  the  i*  instant,  gt{ zj ,  given  by 
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(5-8) 


g,<X)  -  -grW.G)) 


5.1.1  Simplification  of  the  LO  Processor  using  Markov  Processes 

The  LO  processor  of  Eq.  (5-7)  requires  the  N-point  joint  PDF  of  the  baseband  channel 
noise.  In  practice,  it  is  usually  difficult  to  obtain  this  PDF  either  theoretically  or  through 
estimation  techniques.  However,  if  it  is  possible  to  model  the  channel  noise  using  Markov 
random  vectors,  then  generation  of  the  joint  PDF  can  be  greatly  simplified. 


Recall  that  by  using  Bayes  rule,  an  N-point  joint  PDF,  /B(^) ,  of  a  random  vector  ~n  may 
be  written  as  [Woze65]: 


fS  9 fife I  *^2»*  *  *  ’Zn)  ** fifes  I  ^1  »•  *  •  , . . .  ,2V_,) 

~ff-ZN I  I  * ->ZN-l)fife\  ■Zn- 2) 

=/b(ZwI  Z1 » ’ ' '  'ZN-^ffeff-l  I  Zl » ’  •  ’  ’ZN- j)  *  •  ’fifel  I fifed 

However,  if  the  vector  is  an  M*  order  Markov  random  vector,  then 


(5-9) 


ffej\ZV'iZ]-d  ~ffej  I  Zj-M> 


(5-10) 


In  particular,  if  nN  is  a  first  order  Markov  random  vector,  then  ffej\zv...,zj_x)Mfn(zj\zj.l)  and 
Eq.  (5-9)  reduces  to 

(M1) 

j’\ 

with  fn{zx\zj±fn(zx). 
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Using  Eq.  (5-11)  and  assuming  ^  is  a  first  order  Markov  random  vector,  the 
nonlinearity  at  the  i*  instant  in  Eq.  (5-8)  becomes 


g,(Z>  =  -i-mwy) •  -±\n  n/.mv.) 

OZ,  OZf  |_y-l 

az,  j.  i 


(5-12) 


or  equivalently, 

SiO)  =«?Ak-.)  +gfa.t  U<)  ■  (5"13) 

where 

«,(z,iv.)  ^  4ln^iz<-»  (5"14) 

az, 

and 

£,(*,♦.  I z.)  A  I  *<»  (5_15) 

Thus,  for  the  first  order  Markov  case  the  general  N-point  nonlinear  processor,  g((  ,  becomes 

the  simplified  three  point  processor,  ^.(Z/.pZ,,^,),  of  Eq.  (5-13).  Finally,  by  incorporating  the 
first  order  Markov  nonlinearity  of  Eq.  (5-13)  into  Eq.  (5-7),  the  LO  processor  with  memory 
reduces  to  the  following  expression: 
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(5-16) 


*« 


3E  <7 

i*t 


where  ^(z^ZpZj) ^&i(zi,z2)  £^'jv-pzjv*^/v*i)££w(^iv-i,ziv)' 

The  LO  processor  given  in  Eq.  (5-16)  for  the  first  order  Markov  case  may  be  further 
simplified  to  include  only  two  point  nonlinearities,  gf(z,,zM),  by  modifying  the  correlation 

sequence  in  the  following  manner.  Using  Eq.  (5-13)  in  Eq.  (5-16),  L  becomes 


<•1 

Manipulation  of  the  summation  indices  provides  the  following  system  of  equations: 

N  N 

L=E  +E  wMz) 

i*i  <*t 

N  N*  I 

-E  k-.)+E  c,-i  Si-Mi  l  z-«) 


(5-17) 


i*i 

N 


I*  2 


*^nS  Mn*\  I  zn)  ~  C*8<M\  1 2o) 


i*i 


-E 


■•I 


c.  +c.  , 

I  l-l 


(5-18) 


+CffSfM^N*l  I  “  ^0^o(Zl  I  ^o) 


/•  /  V 

Using  Eq.  (5-14)  and  noting  that  /;(z,|z,.l)a-^~~  [Woze65],  ^(zjz,.,)  may  be  written  as: 
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■  -^-(z‘^] 


(5-19) 


or 


SAM  “*&*-.) 


(5-20) 


□ 

where  gi(zitzi_l)&-— In  addition,  a  new  correlating  sequence,  C,  may  be  defined 

3Z;  ^  j 


as 


+c,.  | 


Thus,  the  3-point  test  statistic  of  Eq.  (5-16)  reduces  to  the  following 


2-point  test  statistic  2: 

£<=E  Wzp*i-.)  (5’21) 

(•i 

The  resulting  expression  in  Eq.  (5-21)  is  the  2-point  equivalent  processor  discussed  in 
[Drap89]3. 


2  The  processor  given  by  Eq.  (5-21)  is  derived  from  results  presented  in  [Drap89]  under  the  assumption  that 
cySj,ZN.i  Uw)  =0  “d  C0<?o(Z,  I  Zo)  These  conditions  can  be  met  by  defining  C0  =  0  andg^Z^,  *0 
since  i  =  l N.  This  was  not  explicitly  stated  in  lDrap89J. 

3  There  is  an  error  in  Figure  2-4  of  [Drap89].  The  second  processor  in  this  figure  is  the  equivalent  2-point 
processor.  However,  the  block  diagram  indicates  that  three  points  are  necessary.  Thus,  gi(Zi^vZi,Zi_i)  should  be 
replaced  with  g^Z,.,). 
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5.1.2 


Derivation  of  the  Theoretical  LO  Gain  for  ID  Signalling 


It  is  possible  to  derive  analytical  expressions  for  the  performance  gain  associated  with 
the  LO  processor  with  memory  if,  under  the  central  limit  theorem  [Woze65],  the  statistics  of  L 
are  assumed  to  be  Gaussian.  Then,  the  probability  of  detection,  PD ,  and  the  probability  of  false 
alarm,  PF,  for  the  nonlinear  receiver  are  given  by  [Drap89]: 


VErf 


7-<L,> 


(5-22) 


and 


PF  =  Erf 


y-<L0> 


(5-23) 


where  <L,>  and  aL  are  the  mean  and  standard  deviation,  respectively,  of  L  under  hypothesis 
Hj,  <L0>  and  are  the  mean  and  standard  deviation,  respectively,  of  L  under  H0,  and 

Erf(x)  = — —  [  e  Tdx  [Drap85].  It  can  be  shown  [Drap85,Drap89]  that  Eqs.  (5-22)  and  (5-23) 

are  a  function  of  the  Signal-to-Noise  Ratio  (SNR)  at  the  output  of  the  LO  processor,  SNRTAP  *, 
which  is  given  by 


4  The  term  SNRUP  is  used  here  to  be  consistent  with  [Drap89].  In  [Drap89],  the  LO  processor  with  memory 
is  called  the  Time  Amplitude  Processor  (TAP),  with  corresponding  output  SNR  equal  to  SNRTAP. 
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The  following  is  a  sketch  of  the  derivation  presented  in  [Drap89]  for  the  quantity 
(<£,>-  <L0>).  The  decision  statistic,  L,  is  given  by  5 

where  G  is  used  instead  of  G  to  indicate  that  this  is  the  actual  LO  processor  implementation, 
not  the  theoretically  optimum  one.  The  expression  for  <  L,  >  is  given  by 


5  The  —  scale  factor  is  introduced  in  this  equation  in  [Drap89]  as  a  result  of  implementing  the  required 
N 

summation  for  the  LO  receiver  with  a  standard  linear  correlator.  The  effects  of  this  scale  factor  can  be  removed 
by  appropriately  modifying  the  decision  threshold. 
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(5-27) 


<*■>-}  |  ■Jy hA^S 

z  z 


where yjjW(*)  is  equivalent  to^|f(*)  since  under  hypothesis  Ht  it  is  assumed  that  the  vector cw 
is  transmitted.  Assuming  that  the  signal  is  much  weaker  than  the  channel  interference,  ( • ) 
can  be  approximated  using  a  Taylor  series  expansion  as 


-(1 


(5-28) 


where  Gj^)  is  given  by  Eqs.  (5-7)  and  (5-8).  Substituting  Eq.  (5-28)  into  Eq.  (5-27)  and 
identifying  terms,  it  can  be  shown  that 

<LX>  -  <L0>  |  \cNGf^zt^GN(ztd (5-29) 

ZZ 


or 

<Lt  >  ~  <Lq>  -  i<^Gwr(%>r-  (5-30) 

where  <  •  >  —  indicates  expectation  with  respect  to  ^  (under  Ho)  and  cN.  Next,  by  assuming 
that  the  correlation  matrix  of  ~cN  is  the  identity  matrix  [Drap89],  and  by  using  the  identity 
Py*Tr(7p),  E<1-  (5-30)  becomes 
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(5-31) 


where 


<Lt>  -  <L0>  *lTr(<CvG>-r) 


T^«5Ar>-J 


*5  <u>-‘ 


(5-32) 


The  expression  for  a ^  can  be  derived  as  follows.  Assuming  that  <L0>  »0  because 
‘zN:s~nN  under  H0,  then  »  <Lq>  .  Using  Eq.  (5-26)  and  noting  that  ^  and  are 
statistically  independent  under  H0>  the  expression  for  <L^>  is 

<L0>  *||  Gs  Gn  qj  fn(  zjfc(  cj  d  Zyd  (5-33) 

J,  ^ 

Again  assuming  that  the  correlation  matrix  of  is  the  identity  matrix,  Eq.  (5-33)  reduces  to 

<l°2>  *4Tr<^Ar>T 

(5-34) 

'  N 

where 
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(5-35) 


■it  <«->-, 

w  j»i 


Finally,  by  substituting  Eqs.  (5-32)  and  (5-35)  into  Eq.  (5-25),  the  output  SNR  for  the 
LO  receiver  becomes 


where 


K1 

SNRtap  =s2N~  =s2Nk 

Au 


(5-36) 


(5-37) 


The  expressions  for  Km  and  Kv  of  Eqs.  (5-32)  and  (5-35)  may  be  simplified  if  the  noise 
PDF,  /„('),  is  spherically  symmetric6.  In  this  case,  <^gt>  *  <gjgj>  A  <§g>  and 

<$ii>  3  <gj>  4  <g2>  for  all  i  and  j,  and  Eqs.  (5-32)  and  (5-35)  reduce  to 

Km=<gg>  (5-38) 

and 

K¥ 3  <g2>  (5-39) 


*  In  [Drap89],  the  assumption  of  stationarity,  not  spherical  symmetry,  is  used  to  obtain  the  results  in  Eqs.  (5-38) 

and  (5-39).  However,  stationarity  does  not  guarantee  that  <gtg,>  =  <gjgj>  and  < gf  >  =<|/>  for  all 
i  and  j,  as  will  be  shown  through  simulation  in  Section  5.3. 
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The  gain  of  the  LO  receiver  relative  to  a  linear  receiver  may  be  calculated  by  dividing 
Eq.  (5-36)  by  the  output  SNR  for  a  linear  receiver.  The  output  SNR  for  a  linear  receiver, 

SNRm ,  may  be  computed  by  substituting  ^  for  6^4)  in  Eqs.  (5-32)  and  (5-35),  and  using  the 
result  in  Eq.  (5-36).  This  gives  the  expression 

N 

SNRu  =  LJ1  (5-40) 

On 


where  <?l  is  the  variance  of  the  total'channel  noise.  Dividing  Eq.  (5-36)  by  Eq.  (5-40),  the 
resulting  gain  of  the  LO  receiver  relative  to  a  linear  receiver  is  given  by 


SNRTap  _  2Ki  2 
SNR^  ~°m  Kv  =  °m 


(5-41) 


The  following  is  a  summary  of  the  assumptions  made  in  the  derivation  of  SNRTAP  and  GTAP  : 

•  the  signal  of  interest  is  much  weaker  than  the  total  channel  noise 

•  the  statistics  of  L  are  Gaussian 

•  \ 

•  <L0>  *0,  and  thus  =  <L$> 

•  the  individual  samples  of  cN  are  uncorrelated 

•  the  vectors  4,  and  cv  are  independent  under  hypothesis  H0. 


5.2  Derivation  of  the  LO  Processor  with  Memory  for  2D  Signalling 

Similar  to  the  case  of  ID  signalling,  the  LO  processor  with  memory  for  2D  signalling 
is  derived  in  [Drap89]  as  the  solution  to  the  following  radar  problem:  After  translation  of  the 
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I  and  Q  channels  to  baseband,  consider  the  following  two  hypotheses: 

H{.  xt=Sc os(a+6)+nx 

y.*5  sin(a,+d)  ♦/«„ 

*  (5-42) 

where  (xityi)  is  the  sample  pair  of  the  I  and  Q  channels  at  the  i*  instant,  a,  is  the  i®  sample  of 
the  information  sequence  .(phase  modulation  is  assumed),  S  and  9  are  unknown  amplitude  and 
phase  components  associated  with  the  channel,  and  (nx,ny)  is  the  sample  pair  of  the  total  noise 

in  the  I  and  Q  channels  at  the  i®  instant.  The  following  assumptions  are  made  in  the  derivation 
of  the  LO  processor: 

•  a,  is  independent  of  5  and  9 

•  a(  is  independent  of  nx  and  ny 

•  S  and  9  are  independent  of  nx  and  ny 

•  the  underlying  random  noise  process,  n(t),  which  when  sampled  in  the  I  and  Q  channels 
produces  nx  and  ny ,  is  stationary  and  has  zero  mean 

•  the  signal  of  interest  is  much  weaker  than  the  total  channel  noise. 


The  received  signal  under  hypothesis  H,  of  Eq.  (5-42)  can  be  written  in  a  more  compact  form: 

R^nt  (5-43) 

where  y,]r,  'Si  =  T(ai)l(S,9)  with 


J(S,9)  =  S  e(6)  =S[cos(0)  sin(0)]r 


(5-44) 


representing  the  unknown  parameters,  and 
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I 


tw 


cos  (or.)  -sin(a() 
sinfa,)  cosCa,) 


(5-45) 


representing  the  known  parameters.  The  matrix  T{a)  of  Eq.  (5-45)  has  the  following  useful 
properties: 

•  r'coO-rv,) 

•  n~a)  *  7*r(a() ,  and 

•  nat)n<xJ)*T(c'i+<*j) 


The  LO  processor  is  derived  from  the  standard  likelihood  ratio  test  [VanT68]: 


_  _  Hx 

,  fr\HSRV"^RnW  > 

no 


(5-46) 


where  N  is  the  number  of  samples  per  signal  period.  Using  the  assumptions  outlined  at  the 
beginning  of  Section  5.2,  the  conditional  PDF  /r|H(»)  can  be  written  as 


Jt-7.1  «W|», («!«.)<<« 


4 


(5-47) 


where  fn(>)  is  the  joint  PDF  of  n. ,  i=l,...,N,  and  conditioning  on  H,  implies  that  cr,  is  known. 
Also,  from  Eq.  (5-42),  is  given  by 

/„  =f AJ  (5-48) 

Using  Eqs.  (5-47)  and  (5-48)  in  Eq.  (5-46),  the  likelihood  ratio  becomes 
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f  JS"- '^fy1  _ /  (?i H,)dl 

\  /.(*, . jg 

f '(»/„*,( 


(5-49) 


where 


l(S) a ^ 


(5-50) 


indicates  the  dependence  of  the  likelihood  ratio  on  T,  i.e.  on  the  unknown  amplitude  and  phase 
terms.  Assuming  that  the  desired  signal,  ~S.,  is  much  weaker  than  the  channel  noise,  7£,  the 
expression  for  1(5)  may  be  approximated  via  a  multidimensional  Taylor  series  expansion  as 


/(T)-i*7T-ITa#T 

2 


(5-51) 


where  L  is  a  [2x1]  vector  defined  as 


_  Y.  rr(qj 


(5-52) 


and  M  is  the  [2x2]  matrix  defined  as 


//n  o  \  —  in 


fn(Rv...,RN)  i-i  y-i 


(5-53) 


Since  the  second  order  term,  5,  is  on  the  order  of  (|T|(2  4"<5T5,  M  is  replaced  by  its 


average  value  [Drap89],  M,  given  by 

M=  |  M(Rv...,RN)fll(Rlt...,RN)dRi  ...dRN 
resulting  in  the  following  approximation  of  Eq.  (5-51): 

i(8)  ~  i 

2 


(5-54) 


(5-55) 


5.2.1  Statistics  for  the  LO  Likelihood  Ratio  Assuming  Known  Signal  Amplitude 

and  Phase  Offset 


An  interesting  result  may  be  observed  if  the  likelihood  ratio  of  Eq.  (5-55)  is  examined 
for  a  known  T.  If  T  is  known,  then  the  likelihood  ratio  of  Eq.  (5-49)  reduces  to  1  =  1(5)  and 
the  likelihood  ratio  test  becomes 

1(8)  >7  ^56) 


An  equivalent  test  statistic,  as  discussed  in  [Drap89],  is 


el(1)  <  ey 


(5-57) 


Furthermore,  the  expression  in  Eq.  (5-55)  for  /( 5)  may  be  written  in  a  slightly  different  fashion 
by  "completing  the  square"  to  yield 
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(5-58) 


where,  from  Eq.  (5-60),  fL{L\Hx)  is  a  Gaussian  PDF  with  mean  equal  to  Mb,  fL{L\H^  is  a 
Gaussian  PDF  with  zero  mean,  and  both  have  a  covariance  matrix  equal  to  M.  In  this  context, 

the  decision  statistic  L  may  be  considered  to  be  a  Gaussian  random  vector  under  both  the 
hypotheses,  H,  and  H0.  Thus,  the  development  of  Eq.  (5-60)  provides  another  argument,  in 
addition  to  the  central  limit  theorem,  that  the  statistics  of  Z  are  Gaussian. 
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5.2.2 


Development  of  the  LO  Likelihood  Ratio  Assuming  an  Unknown  Nonrandom 
Signal  Amplitude  and  a  Random  Phase  Offset 


The  most  useful  scenario  to  examine  is  when  the  signal  amplitude  is  unknown,  but  not 
necessarily  random,  and  the  phase  offset  is  a  random  variable.  In  this  case,  H,  becomes  a 
composite  hypothesis  [VanT68]  because  the  overall  problem  is  still  a  binary  detection  problem, 
but  the  received  signal  under  H,  is  a  function  of  a  random  variable,  6 ,  and  an  unknown  albeit 
nonrandom  parameter  5.  Assuming  S  is  a  known  parameter  for  the  moment,  the  likelihood  ratio 
for  this  scenario  is  given  by  [VanT68] 


fr\ W,(^l » * "  • » ^0  I  »o) 

_  [/„„,(*, . SylMf, )fmV>  I 

/„*(*«  *  *  •  •  »^/V  I  tfo) 


(5-62) 


Now,  for  the  case  when  S  is  unknown,  it  is  intuitively  obvious  that  the  detector  with  the  best 
performance  will  be  one  that  is  able  to  measure  S  perfectly  and  then  use  this  value  in  the 
optimum  likelihood  ratio  test.  However,  in  most  cases  this  type  of  detector  is  unrealizable  due 
to  errors  caused  by  the  channel  noise.  Thus,  one  logical  method  for  constructing  the  likelihood 
ratio  is  to  first  errimate  S  and  then  use  this  estimate  in  the  likelihood  ratio  test  as  if  it  were 
exact.  One  possible  estimator  for  S  is  the  Maximum  Likelihood  Estimate  (MLE),  denoted  by 

5^.  The  value  of  Sml  is  equal  to  the  value  of  5  which  maximizes /r|i(/?,,...,/?Ar|5) 7,  or  in  other 

words,  is  the  value  of  S  most  likely  to  have  produced  the  given  value  of  the  received  signal 

vector.  Using  the  above  argument,  one  can  construct  the  Generalized  Likelihood  Ratio  (GLR) 
[VanT68],  which  in  this  case  is 

7  This  is  the  reason  why  one  of  the  main  focal  points  of  Section  3. 2.2.2. 1  of  [Drap89]  is  determining  the  value 
of  S  which  maximizes  the  decision  statistic,  /. 
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(5-63) 


^  I  J-S^ 

,a  /^h(R1,...,Rn\H<) 


maX/r|S.//,(^..-^l^|) 

iS 

X)  h(^1  *  •  •  •  I  *0> 


Using  the  development  of  Eq.  (5-63),  Eq.  (5-62)  can  be  modified  to  yield  the  GLR  for  the  case 
of  unknown  nonrandom  signal  amplitude  and  random  phase  offset 


|X|«  S  tf  !*•••» 

fr\H^\ » *  *  *  *  I  "o> 


(5-64) 


Also,  /mH(  •)=/«,„(•)  since  S  and  d  are  assumed  to  be  statistically  independent.  Thus 


/..max  | ^ 


.max 

O  I 


(5-65) 


=  <  max  l(S,9)>e 

s 


where  <  •  > ,  indicates  expectation  with  respeu  to  6 , 


Substituting  Eq.  (5-44)  into  Eq.  (5-51)  yields  the  following  LO  expression  for  1(5): 
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(5-66) 


1(d)  =/(S,0)  =  1  +s7(d)rL-is2e(d)TM7(6) 


Provided  that  M  >  0 ,  the  maximum  of  l(S,d)  with  respect  to  S  occurs  when 


5=5‘(g)a-^-?  - 

e(6)TMe(6 ) 


(5-67) 


Substituting  Eq.  (5-67)  into  Eq.  (5-66)  produces  the  expression 

max  l(S,d)  =l(Sm,0)*l  +  . 

5  2  e(d)TMe(d) 

Using  Eq.  (5-68)  in  Eq.  (5-65)  yields  the  generalized  likelihood  ratio  test 


(5-68) 


<i(e(6)TL)2>,  >  y<e(6)TMe(e)>, 


(5  -69) 


where  7=7-1.  Using  the  identity 


e(9)TM e(d) > s=1x[m  <  e(0)  e(9)T>  J 


(5-70) 


and  assuming  that  9  is  uniformly  distributed  over  the  interval  [0,2v)  yields  the  result 


e(9)TM  e(9)  > ,  =Tr[M  <  M  W)T>  J 


(5-71) 


since  <  e(9)  e(9)T>9  is  a  two  dimensional  identity  matrix.  Similarly,  it  can  be  shown  that 


<(e(0)TL)2>i  =  LTL&  |iJ2 


(5-72) 


Substituting  Eqs.  (5-71)  and  (5-72)  in  Eq.  (5-69)  gives  the  expression 
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(5-73) 


_  W> 

\ 

n0 

Finally,  recalling  the  discussion  in  Section  5.2.1,  M  may  be  though  of  as  the  covariance  matrix 
of  L ,  and  as  such 

Tr{iWr]*Tr<ZTr>  j 

= <rr>*  (5-74) 

A<||I0l2>3f 

where  <  •  indicates  expectation  with  respect  to  the  received  signal  under  Hq  (the  channel 
noise)  *.  Thus,  the  resulting  LO  likelihood  ratio  test  for  the  case  of  an  unknown  nonrandom 
signal  amplitude  and  a  random  phase  offset  is 

1  INI2  >y  (5-75) 

2<\\L0\\2>  * 


5.2.3  Derivation  of  the  Theoretical  LO  Gain  for  2D  Signalling 

The  performance  measures  of  interest,  as  in  the  case  for  ID  signalling,  are  PD  and  PF. 
Similar  to  the  ID  signalling  case,  and  discussed  in  [Drap85],  the  output  signal  to  noise  ratio, 

*  This  result  is  presented  in  [Drap89].  Note  that  for  Eq.  (5-74)  to  be  true,  the  covariance  matrix  of  L  under 

H,  must  be  approximately  equal  to  the  covariance  matrix  of  L  under  Ho,  and  the  mean  of  L  under  H*  must  be 
approximately  equal  to  the  zero  vector.  These  assumptions  were  not  explicitly  stated  in  [Drap89]. 


SNRtap ,  of  the  LO  processor  determines  the  performance.  SNRT4J,  is  given  by  the  expression 
[Drap893’ 


snrtap*e\ 


<L,j  5>  -<Lq\  8>  112 

<il!J2> 


(5-76) 


where  <  •  >  represents  expectation  with  respect  to  the  information  sequence  and  the  received 
signal  under  hypothesis  H0  or  H,  as  required,  and  £{  •  }y  indicates  expectation  with  respect  to 


5. 

As  indicated  in  Eq.  (5-76),  the  terms  required  to  evaluate  SNRTAP  are  calculated  given  T, 
and  then  the  result  is  averaged  with  respect  to  S,  or  0 ,  or  both  depending  on  the  given  scenario. 
Recalling  Eq.  (5-52),  the  expression  for  L  may  be  written  as 


T\ct)Gi(Rv...JN) 


(5-77) 


where  Gt(Rv...,RN)  represents  the  implemented/approximated  version  of 


n(~Q  1?  \  ^  _  1  d/,(£,,...,/?v) 

■  — *  — 

•  *  '  /•  /  n  n 


(5-78) 


Rn)  dRi 


in  Eq.  (5-52).  Since  ait  5,  and  nt  are  independent,  and  since  the  received  signal  under 

hypothesis  H,  is  a  function  of  all  three  of  these  variables,  the  expression  for  <  L  |  5  >  may  be 
written  as 


9  Note  that  in  [Drap89]  the  required  conditioning  on  5  was  omitted. 
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(5-79) 


<L||  5>  |  LfKa(Rl,...,RN,aN\  5 )dRr..dRNdaN 

Lf*(R\  •  >  *v  -  SN)fa(  aN)dRl . .  .dRydotjy 


<\  R 


a 


where  aN  =  [afl,...,ajV]7'.  Assuming  that  the  signal  is  much  weaker  than  the  channel  noise,  a  first 
order  multidimensional  Taylor  series  expansion  of  fn(  • )  may  be  used,  resulting  in  the  expression 

. M  -/.(A,,...,**) 

where 


GjiRv-'RjnccjUie) 


(5-80) 


(5-81) 


Using  Eqs.  (5-77)  and  (5-80)  in  Eq.  (5-79),  and  identifying  the  term  that  represents 
<L0|  5>  »  <L0>  ,  yields  the  following  result 


<L,|5>  -  <JL0|5>  ~SNKm  e(6)  (5-82) 

where  Km  is  the  [2x2]  matrix  defined  as 


Ar-*iEE<?’rWC,5/7to!))>,- 

i«i  /» i 


(5-83) 


where  <  •  >*-  indicates  expectation  with  respect  to  the  received  signal  and  as.  Since  in  the 
approximation  of  Eq.  (5-80),  the  received  signal  is  assumed  to  be  nearly  equal  to  the  noise,  and/t 
and  a.  are  assumed  to  be  independent,  Eq.  (5-83)  may  be  written  as 
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(5-84) 


K  -  4E  t  <  TT(a)  <  <5,  Gj>  tlX*)  >  - 

W  i»l  y»l 

Finally,  if  it  is  assumed  that  a(.  is  independent  of  a;,  and  that  a,  is  uniformly  distributed  over 
the  region  (0,2ir),  then  the  expression  for  Km  simplifies  to 


k.-^E  <rr(a,)  <G1C,r>, 7tc,)>- 
n  «•! 


(5-85) 


The  expression  for  <  L0 1| 2  >  is  given  by: 


lirj2>  =  <  II ll 2 1 s>  -  \ }  (5_8(S) 


1 


Substituting  Eq.  (5-77)  into  Eq.  (5-86)  yields  the  following  expression  for  <  ||L0||2> 


<  l0\2>=nkv 


(5-87) 


where  Kv  is  the  scalar  parameter  defined  as 


N  <V 


*uiE£  <clna)T\ai)dl>^ 

I* i  >*i 


(5-88) 


Again,  since  n  is  independent  of  a.,  and  if  it  is  assumed  that  a(.  is  independent  of  a;,  with  a. 
uniformly  distributed  over  the  region  [0,2ir) ,  the  Eq.  (5-88)  reduces  to 


*.-^E  <G,rGl>i 


(5-89) 


The  expressions  for  Km  and  Kv  of  Eqs.  (5-85)  and  (5-89)  may  be  simplified  if  the  noise 


PDF,  /,(•)  is  spherically  symmetric  10.  In  this  case,  <GfiJ>  *  <GjGj>  a  <0Gt>  and 
<d^dt>  »  <(jJr6j>  a  <GtC>  for  all  i  and  j,  and  Eqs.  (5-85)  and  (5-89)  reduce  to 

K  =  <GGt>  (5-90) 

n 

and 

A><GrG>  (5-91) 

% 

Finally,  substituting  Eqs.  (5-82)  and  (5-87)  into  Eq.  (5-76)  yields  the  following 
expression  for  SNRTAP : 

SNRtap  .  £e{s'  II  KjH)  II  ’}T  (5-92) 

For  the  case  when  5  is  deterministic  and  6  is  uniform  over  [0,2t)  ,  Eq.  (5-92)  becomes 

or 

=  r  KmK;±.l7?-de  ,5‘94) 

However,  since  0  is  uniform  over  [0,2tt) 

10  In  [Drap89],  the  assumption  of  stationarity,  not  spherical  symmetry,  is  used  to  obtain  the  results  in  Eqs. 
(5-90)  and  (5-91).  However,  stationarity  does  not  guarantee  that  <  G(  Gi  >  =  <  Gj  G,  >  and 
<  G(  G;>  *  <  G;  G;>  for  all  /  and  j,  as  will  be  shown  through  simulation  in  Section  5.3. 
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Thus,  Eq.  (5-94)  reduces  to 


M 


2i 

f  7~e  Td8  =  ]- 

1  0 

i  2 

0  1_ 

T'fcX 


smTAr*s'N  2K 


(5-95) 


(5-96) 


The  theoretical  gain  of  the  LO  processor  relative  to  a  linear  receiver,  GTAf,  may  be 
computed  by  taking  the  ratio  of  SNRTAP  to  the  output  SNR  of  a  linear  receiver,  SNR^.  The 

expression  for  SNR^  is  calculated  by  using  Gi  =  [jc(  y.]r  to  evaluate  Km  and  Kv  in  Eq.  (5-96). 
Evaluating  <d,G/>y  in  Eq.  (5-84)  for  Gi  =  [Jt1,y1]r  yields 


<  GfiJ>j  =  -  |  Q(R} RJ  dR, dRN 

I 


(5-97) 


where 


Q(  Rn) 


„  dfn{Rv...,Rs)  x  dfn(Rlt...,RN) 


dx, 


dyj 


yt 


dfH(Rit...,RN)  dfn(R\ . Rn) 


dXj 


yt- 


dyj 


(5-98) 


Evaluating  the  two  off-diagonal  terms,  and  assuming  that  the  tails  of  /„(  • )  approach  zero  as  x} 
and  yj  approach  ±  oo ,  yields 
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I 


y  =oo  __ 

'  dRr.dR  dXjdR  r.dRw= 0 

y  =-oo 

«**  Ji  '  (5-99) 


and 


dfn(Rv...,RJd-^  d-  s_ 


dyj 


i 


_  _  -  *7*°°  _  ^ 


except  x. 


yt  fSRx . R*)\  dRr.dRJ.ldyjdRj.l..dRfi= 0 

(5-100) 


*r-°° 


For  i-j,  using  integration  by  parts  the  diagonal  terms  of  Eq.  (5-97)  reduce  to 

_  □ 


Qf  _  _  _  X^*oo 

-\x~dRv..dRs  =  -  f  xtfm(Rv...M¥  dyr.dyMdyi.r.dyN+(fn(Rv..RJdRr.. 


dR K 


except  yt 


,-..,RN)dRv..dRs 


=  1 


(5-101) 


and 
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again  since  the  tails  of  the  channel  noise  PDF  are  assumed  to  decay  to  zero.  Substituting  the 
results  of  Eqs.  (5-99)  through  (5-104)  into  Eq.  (5-84),  and  noting  that  TT(a)I\a)  is  the  two 
dimensional  identity  matrix,  K  becomes 
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(5-105) 


where 


1  0 


(5-106) 


To  evaluate  Kv^,  G(  =  [*1,y.]r  is  substituted  into  Eq.  (5-88)  which  yields  the  expression 

K^aUal-al  (5-107) 

where 

<£  =  | xffH(Rx ,...,RS) dRl . . .dRs  , 

°)  3  |  ytfn&i  »••• to  <&X -dRy  , 

Substituting  Eq.  (5-105)  and  (5-108)  into  Eq.  (5-96)  yields  the  following  result  for  SNRm  : 

(5-109) 

On 


for  all  / 

for  all  i 


(5-108) 


Finally,  dividing  Eq.  (5-96)  by  (5-109)  yields  the  expression  for  GTAP  : 


(5-110) 
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5.2.4  Theoretical  LO  Output  SNR  and  Gain  for  BPSK  Modulation 

The  derivation  of  SNRTAP  and  GTA?  for  the  case  of  BPSK  modulation  is  examined  in  this 
section.  For  BPSK,  a(.  ®0°  or  180°  and  7(a)  is 


where 


Also,  Kv  becomes 


where 


Ks\  G&f&yJd&y* 

Kmm  =  {  GyGy, 

=  {  ^G,  fn(Zy,Jjdx^dyN 

*  |  GxGy,  dXydyy 

J  Gxyn(vN)d%AyN 
\m\  GyXiWj&Jy* 


(5-11(5) 


(5-117) 


(5-118) 


Substituting  Eqs.  (5-115)  and  (5-117)  into  Eq.  (5-96)  yields  the  resulting  expression  for  SNRTAP 
for  BPSK: 


s2  EE^  K  +£ E K*  K~  +  E E K*  Km  +E E K~  K« 

•Vf-,  4?  4^  %  4?  "  m*.  %  4f  4*  "S..  m* 

_  i»l  7»l  i»l  ;»l  i»l  ./*!  i*I  /*1 

%  <v*4 


(5-119) 


Similarly,  Gr<4/,  for  the  case  of  BPSK  modulation  is: 
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If  the  channel  noise  PDF,  /„(•),  is  spherically  symmetric,  the  expressions  for  SNRtaP  and 
GTap  may  be  greatly  simplified.  In  this  case,  Km  =Km  &Km,K  =Km  &Km  ,  Km  =Km  4  AT  , 

“i  m  mnt  mn  moi  mn 

Km  -Km  &K  ,  and  also  K  -Kv  4  K  and  Kv  =KV  for  all  /  and  j.  Thus,  Eqs.  (5-119) 


and  (5-120)  simplify  to 


SNRtap  = 


2(*.>M 


(5-121) 


(5-122) 


5.3  Simulation  Results  for  the  LO  Processor  with  Memory 

To  provide  insight  into  the  RDAT  theory  and  methodology,  and  to  provide  a  means  of 
verification,  an  interference  scenario  discussed  in  [Drap89]  was  simulated  using  the  MATLAB™ 
software  package  by  The  Math  Works,  Inc.  The  scenario,  depicted  in  Fig.  (5-1),  involves  ID 
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signalling  and  a  channel  interference  consisting  of:  i)  additive  white  Gaussian  noise,  w(n),  with 
variance  given  by  ol,  and  ii)  colored  Gaussian  noise,  j(n),  with  variance  given  by  a].  The 
sequence  j(rt)  is  obtained  as  the  output  of  an  N-point  Moving  Average  (MA)  process,  whose 
input  is  white  Gaussian  noise,  v(n),  having  variance  equal  to  ac.  Thus,  the  colored  Gaussian 
noise  sequence,  j(n),  is  given  by  the  equation 

j(n)  =  -L  53  v(n  -k)  (5-123) 


The  autocorrelation  sequence,  \(k)  =E{n(n)n(n+k)}  =  <  n(n)n(n+k)  >  ,  of  the  total  channel  noise, 
n(n)  =j(n)  +  w(n) ,  is  given  by 


ol  +  Oc  . 
N-\k\ 
N 

0  , 


2 

Oc  > 


l£|*|  <N 
k^N 


(5-124) 


Finally,  the  covariance  matrix  for  the  total  channel  noise  is  given  by 


'  X(0) 

X(l) 

X(2) 

-  \(P~  1) 

X(1) 

X(0) 

X(l) 

•••  \(p-  2) 

A  = 

X(2) 

X(l) 

X(0) 

-  \(p- 3) 

X(p-l) 

\(p-2) 

X(/7-3) 

-  X(0) 

(5-125) 


where  p  is  the  number  of  received  samples  processed  by  the  LO  processor  per  signal  period. 
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LO  Processing  with  Memory 
Colored  Gaussian  Noise  Scenario 


Figure  (5-1) 


Since  the  total  channel  noise  is  the  sum  of  two  zero  mean  Gaussian  noise  processes,  its 
PDF  is  also  Gaussian  and  is  given  by 


/fl(")=(2xy’/2|A|-,/2e  7 


(5-126) 


where  n=[nx  n2  ...  np]T.  Thus,  the  nonlinearity  at  the  i*  instant,  g„  is  given  by 


gr-j^Wn(n))  =  nrA-'  et 


(5-127) 


where  ~et  is  a  vector  of  length  p  whose  elements  are  all  zero  except  for  the  i*  element  which  is 
equal  to  one.  Recalling  that  for  ID  signalling  the  gain  of  the  LO  processor  relative  to  a  linear 
receiver,  GTAP ,  is  given  by 
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(5-128) 


^tap  =  kan 
K* 

where  A:  =  — 2 ,  Km  and  K,  are  given  by  Eqs.  (5-32)  and  (5-35),  a\  =  a\*aw  in  this  scenario,  and 

that  for  the  LO  receiver  gt  =gr  Thus,  the  value  of  k  for  the  LO  receiver  may  be  computed  as 

*  =  (5-129) 

ptl 

Substituting  Eq.  (5-127)  into  Eq.  (5-129),  and  noting  that  <H~nT>  =A,  yields  the  following 
result  for  k: 


P  i-i 


Finally,  the  expression  for  GTap  is  given  by 


«»-7 

P  i»l 


(5-130) 


(5-131) 


Figures  (5-2)  to  (5-5)  illustrate  the  resulting  values  of  GTAP  in  relation  to  the  LO 
processor  length,  p ,  and  the  length  of  the  colored  Gaussian  MA  process,  N,  for  ac-a\-  1, 

a)  -  10a«  =  10 ,  a]  =  100a^,  =  100 ,  and  a]  =  lOOOcrJ,  =  1000 ,  respectively.  As  may  be  seen,  the  gain 

of  the  LO  processor  increases  as  its  length  increases.  It  is  shown  in  [Drap89]  that  as  the  LO 
processor  length  increases,  the  gain  relative  to  a  standard  linear  receiver  approaches  that  of  the 
optimal  filter  for  this  scenario,  the  Wiener  filter.  Furthermore,  for  a  given  LO  processor  length, 
the  gain  increases  as  the  value  of  N  increases.  This  indicates  that  the  LO  processor  is  better  able 
to  remove  the  effects  of  the  colored  Gaussian  noise  as  adjacent  samples  become  more  correlated. 
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In  comparing  the  results  shown  in  Figs.  (5-2)  through  (5-5)  to  those  in  [Drap89]  (see 
Figs.  2-8  to  2-1 1  on  pp.  2-33  to  2-36),  it  was  noticed  that  the  results  were  not  identical.  The 
reason  for  this  discrepancy  is  that  the  results  in  [Drap89]  were  generated  using  the  assumption 
that  since  the  channel  noise  process  was  stationary,  the  value  for  k  could  be  calculated  using  just 
one  gi,  i.e.  using  the  equation 

k-<g?>  (5-132) 

for  any  /.  However,  as  discussed  in  Sections  5.1.2  and  5.2.3,  stationarity  is  not  a  sufficient 

condition  to  ensure  that  <  g,2  >  =  <g/>  for  all  i  and  j.  Instead,  the  PDF  of  the  channel  noise 

must  be  spherically  symmetric,  which  is  not  the  case  in  this  scenario.  To  show  that  stationarity 
is  not  a  sufficient  condition  to  give  the  result  in  Eq.  (5-132),  the  values  of  GTAP  in  relation  to  p 

and  N  were  plotted  for  the  case  of  a \  =  100a?„  =  100  using  Eq.  (5-132)  with  /  =  /.  The  results  of 

this  experiment  are  shown  in  Fig.  (5-6)  and  are  almost  identical  to  the  results  presented  in 
[Drap89]  11 .  It  may  be  observed  that  Figs.  (5-6)  and  (5-4)  are  not  the  same,  but  should  be  if 
stationarity  is  sufficient  to  ensure  the  validity  of  Eq.  (5-132).  Thus,  stationarity  is  not  a 
sufficient  condition  to  warrant  the  simplification  of  the  expression  for  k  given  in  Eq.  (5-129)  to 
that  given  in  Eq.  (5-132). 


11  There  are  some  slight  differences  in  the  plot  shown  in  Fig.  (5-6)  and  that  of  Fig.  2-10  of  [Drap89],  These 
differences  are  most  likely  a  result  of  differences  in  the  resolution  of  the  computer  software  used  to  generate  and 
plot  the  results. 
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Colored  Goussion  N©i,.  Memory  Length  n 
°c  =  <TW2  a  1 

Figure  (5-2) 
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Figure  (5-3) 


LO  Processor  Gain  for  the  Colored  Goussion  Noise  Scenario 


Figure  (5-4) 


LO  Processor  Gain  for  the  Colored  Gaussian  Noise  Scenario 


Figure  (5-5) 
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t0  Processor  Coin  for  the  Colored  Goussion  Noise  Scenono 


Figure  (5-6) 
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6. 


FUTURE  OBJECTIVES 


While  the  work  presented  herein  provides  an  assortment  of  clues  into  the  fundamental 
behavior  of  the  LO  processor,  there  is  still  much  research  to  be  done  in  this  area.  Previous 
sections  of  this  report  addressed  the  memoryless  LO  processor  and  the  many  algorithms  used 
to  implement  it.  Software  simulation  is  the  method  of  choice  for  characterizing  the  performance 
of  this  processor.  Two  simulations  were  presented  which  provided  a  means  of  investigating  the 
Pb  of  the  memoryless  LO  processor  in  a  variety  of  jamming  environments.  The  LO  processor 
with  memory  was  also  discussed,  with  an  emphasis  placed  on  its  derivation  and  the  development 
of  the  theoretical  expressions  for  the  output  SNR  and  gain  relative  to  a  linear  receiver.  A 
number  of  logical  extensions  of  the  research  presented  in  this  report  may  be  identified.  Some 
of  these  include: 

•  Simulation  of  the  LO  processor  with  memory  and  its  incorporation  into  the  QPSK  DSSS 
system.  By  comparing  the  results  of  this  simulation  to  the  memoryless  LO  processor 
results,  the  performance  enhancement  gained  through  the  use  of  memory  can  be 
quantified. 

•  Investigation  of  other  modulation  techniques,  particularly  Continuous  Phase  Frequency 
Shift  Keying  (CPFSK)  [Leib93],  a  form  of  Minimum  Shift  Keying  (MSK),  and  their 
effects  on  the  performance  of  the  LO  processor. 

•  Simulation  of  the  parametric  implementation  of  the  FSA  algorithm.  All  of  the  LO 
processor  implementation  methods  mentioned  in  this  report  employ  a  histogram 
approximation  of  the  received  signal  PDF  at  some  point  in  their  respective  algorithms. 
Since  it  was  seen  that  all  the  methods  had  roughly  the  same  performance  results,  one 
may  conclude  that  the  histogram  is  the  limiting  factor.  By  examining  the  parametric 
implementation  of  the  FSA  algorithm,  it  may  be  possible  to  determine  if  this  hypothesis 
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is  true,  or  if  the  method  of  LO  processor  implementation  has  relatively  little  effect  on 
its  overall  performance. 

•  Adaptation  of  the  Continuous  Polynomial  Approximation  (CPA)  algorithm  [Grim93]  for 
use  as  an  LO  processor  implementation  method.  Through  simulation  the  performance 
of  this  algorithm  can  be  compared  to  the  other  methods  (Histogram,  EBH,  FSA,  and 
MIPA)  previously  discussed. 

LO  processing  is  a  rich  area  of  study.  Intuitively,  the  LO  processor  with  memory  should 
be  more  robust  than  the  memoryless  processor  and  thus  able  to  mitigate  a  wide  array  of  different 
jammer  types,  including  those  for  which  the  memoryless  LO  processor  performed  poorly.  Also, 
whether  or  not  different  LO  processor  implementation  algorithms  have  a  strong  influence  on 
overall  performance  remains  to  be  seen.  The  work  presented  in  this  report,  along  with  the 
future  tasks  mentioned  in  this  section,  will  lead  to  the  definition  of  a  communications  system  that 
is  robust  and  able  to  provide  a  secure  and  reliable  information  exchange  in  a  multitude  of 
interference  scenarios. 
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